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ABSTRACT 


A  mean  streamline  analysis  of  the  effect  of  stator  blade  orienta- 
tion on  the  performance  of  an  axial  flow  compressor  was  performed  by 
means  of  a  computer  program.   Measurements  were  made  on  a  3-stage  axial 
flow  compressor  at  the  Naval  Postgraduate  School  at  six  stator  stagger 
angles  between  23.8  and  44.3   for  a  fixed  orientation  of  the  rotor 
blades.   Maximum  efficiency  and  pressure  ratio  were  measured  at  a 
stator  stagger  angle  of  31.8  .   Results  at  other  blade  settings  showed 
that  by  varying  stator  stagger  angle  with  flow  rate  optimum  efficien- 
cies and  pressure  ratios  can  be  achieved  over  a  wide  range  of  operating 
conditions. 

The  results  of  the  analysis  were  compared  with  the  measured  results, 
Suggestions  are  made  for  improving  the  manner  of  adapting  cascade  test 
data  to  performance  predictions. 

By  applying  a  non-dimensional  deflection  coefficient  it  could  be 
shown  that  minimum  work  input  corresponded  to  maximum  efficiency. 

The  test  compressor  has  a  tip  diameter  of  36  in.  and  a  hub/tip 
ratio  of  0.6.  The  blading  tested  is  of  the  free-vortex  type  with  a 
design  degree  of  reaction  of  0.5.   Tip  speed  was  about  185  ft/sec. 
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SECTION  1 
INTRODUCTION 

Axial  flow  compressors  have  a  wide  range  of  application.   They 
are  usually  matched  to  a  turbine  in  a  set.   Both  open-cycle  sets, 
such  as  aircraft  turbo-jet  engines,  and  closed-cycle  sets  for  power 
generation  are  in  service.   In  every  application  compressor  efficiency 
is  the  critical  factor  in  set  performance.   Accurate  predictions  of 
off-design  compressor  performance  must  be  made  during  the  design  stage 
to  be  able  to  make  corrections  prior  to  manufacturing. 

Theoretical  analysis  of  the  flow  in  an  axial  compressor  involves 
three-dimensional  partial  differential  equations  of  a  complex  nature  (1). 
In  these  equations  it  is  difficult  to  account  for  real  gas  effects  such 
as  boundary  layer  growth  on  machine  walls  and  blade  surfaces.   The 
resulting  wakes  behind  blade  rows  create  non-uniform  conditions;  hence, 
it  becomes  necessary  to  base  prediction  methods  on  experimental  data. 

These  data  would  be  obtained  best  on  an  actual  compressor.   How- 
ever, most  compressors  have  relatively  short  blades  and  small  flow 
annuli.   Even  if  only  small  pressure  probes  were  inserted  between  the 
rows  of  blades,  the  flow  in  these  machines  would  be  substantially 
altered.   Test  rigs  must  therefore  be  used  that  have  large  dimensions. 

A  first  approximation  of  the  flow  in  axial  turbomachines  can  be 
obtained  in  rectilinear  cascade  test  rigs.   The  intersections  of  a 
co-axial  cylinder  with  the  blades  of  a  row  produce  a  series  of  iden- 
tical and  identically  oriented  profiles  which  are  unwrapped  into  a 
plane  to  establish  the  corresponding  rectilinear  cascade.   A  finite 
number  of  straight  blades  with  profiles  similar  to,  but  larger  than 
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those  obtained, is  then  arranged  in  a  rectangular  duct  through  which 
air  is  blown  at  the  appropriate  inflow  angle.   Flow  surveys  taken  ahead 
of  and  behind  the  blade  row  give  experimental  data  of  air  turning  angle 
and  total  pressure  loss  for  various  blade  geometries.   Such  a  recti- 
linear cascade  is  in  use  at  the  Naval  Postgraduate  School.   The  results 
of  extensive  testing  of  compressor  airfoil  shapes  at  NASA  have  been 
summarized  by  Lieblein  (2). 

The  usual  problem  in  compressor  design  consists  in  selecting  flow 
areas,  blade  shapes  and  blade  layout  to  satisfy  design  specifications. 
Additional  studies  are  frequently  necessary  if  the  design  performance 
is  not  reached.   An  example  is  the  study  of  Vavra  which  was  made  when 
the  so-called  Clark  CSN-1  compressor  of  the  ML-1  nuclear  gas  turbine 
of  Aerojet-General  Nucleonics  failed  to  perform  as  required.   Vavra 
analyzed  the  design  changes  proposed  by  the  manufacturers  and  made 
further  recommendations  (3).   However  the  ML-1  project  was  cancelled 
before  the  improved  compressor  could  be  built  so  that  it  was  not 
possible  to  verify  the  suggested  changes  by  experiments. 

The  so-called  inverse  problem  consists  in  predicting  the  perform- 
ance of  an  existing  machine  with  a  known  blading  by  means  of  rectilinear 
cascade  data  and  other  experience  factors.   Such  an  approach  was  carried 
out  by  Gibbons  and  Bartels  (4)  for  the  12-stage  Allis-Chalmers  axial 
flow  compressor  which  supplies  air  to  the  turbine  test  facilities  at 
the  Naval  Postgraduate  School.   Because  of  the  unorthodox  performance 
of  the  first  stages,  and  also  due  to  small  flow  annuli  and  short  blade 
heights,  it  was  not  possible  to  predict  the  performance  of  the  machine 
accurately.   However,  the  analysis  gave  important  indications  for 
design  changes  to  improve  the  performance  of  the  compressor. 
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This  thesis  makes  another  attempt  to  solve  the  inverse  problem 
for  a  3-stage  axial  flow  compressor  built  by  the  California  Institute 
of  Technology  with  the  support  of  the  Office  of  Naval  Research.   The 
compressor  is  sufficiently  large  so  that  the  insertion  of  pressure 
probes  causes  relatively  small  flow  perturbations.   The  machine  was 
designed  to  permit  variations  in  blade  shapes,  blade  angles,  tip 
clearances,  and  staging.   The  installation  at  the  California  Institute 
of  Technology  is  described  by  Bowen,  e_t  a_l.  (5)   This  report  also 
describes  a  theory  of  perfect  fluid  flow  in  axial  flow  turbomachines. 
The  results  of  this  theory  are  compared  with  the  measured  data  of  the 
first  stage  of  a  so-called  "free-vortex"  type  blading.   Part  2  of  the 
report  by  Bowen,  e_t  a_l.  (6),  gave  a  detailed  investigation  of  multi- 
stage flow  for  both  "free-vortex"  blading  and  a  more  complex  "solid 
body  rotation"  blading.   Measurements  of  blade  skin  friction  losses 
were  greater  than  expected  from  cascade  tests.   The  growth  of  boundary 
layer  along  the  inner  and  the  outer  annulus  walls  of  the  flow  was  less 
than  expected.   Alsworth  and  Iura  (7)  conducted  extensive  and  detailed 
measurements  of  flow  patterns  in  a  single  stage  of  "free-vortex" 
blading.   They  carried  out  accurate  measurements  of  the  blade  skin 
friction  losses  and  the  radial  distribution  of  work  input.   The  final 
report  from  the  California  Institute  of  Technology  was  a  hot-wire 
anemometer  study  of  compressor  stall  by  Iura  and  Rannie  (8).   They 
observed  that  stalled  flow  regions  rotate  in  the  direction  of  blade 
rotation  without  changing  shape  but  with  a  speed  that  is  not  propor- 
tional to  rotor  speed. 

After  the  Turbo- Propulsion  Laboratory  of  the  Naval  Postgraduate 
School  was  built,  the  compressor  was  relocated  there.   Since  then 
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it  has  been  used  for  laboratory  courses  to  supplement  instruction 
in  basic  theories  of  turbomachines. 

The  present  study  was  conducted  for  this  compressor  because  of  the 
possibility  of  changing  its  blade  angles.   For  fixed  stator  and  rotor 
blade  angles  a  highly  peaked  curve  of  efficiency  versus  flow  rate  is 
usually  obtained  with  axial  flow  compressors.   At  a  slightly  changed 
stator  blade  angle  the  peak  will  be  displaced  somewhat.   A  family  of 
efficiency  curves  for  various  stator  stagger  angles  is  expected  to 
yield  a  flatter  efficiency  curve  over  a  wider  range  of  flow  rates  for 
fixed  rotor  blade  angles.   The  wide  variation  in  flow  rates  required 
of  the  engines  for  the  supersonic  transport  has  led  to  design  pro- 
posals by  General  Electric  for  compressors  where  the  stator  blade 
angles  can  be  changed  during  operation  of  the  jet  engine.   Hence,  the 
thesis   topic  is  of  current  interest. 

For  the  tests  it  was  necessary  to  obtain  extremely  accurate 
pressure  measurements.   An  extremely  sensitive  bourdon  tube  pressure 
indicator  was  added  to  the  instrumentation.   A  permanent  Pitot-static 
tube  was  installed,  and  accurate  flow  rate  calibrations  were  carried 
out.   A  parallel  theoretical  effort  produced  predictions  of  com- 
pressor performance  by  an  existing  computer  program  which  was  modified 
and  adapted  for  use  on  the  I.  B.  M.  Model  360  computer  system  installed 
at  the  Naval  Postgraduate  School. 

The  author  wishes  to  express  a  debt  of  gratitude  to  the  faculty 
of  the  Department  of  Aeronautics  of  the  Naval  Postgraduate  School  for 
his  engineering  education.   Particular  thanks  are  due  to  Professor 
M.  H.  Vavra  for  his  enthusiasm  and  guidance  in  showing  what  can  be 
done  with  that  education. 
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SECTION  2 
0.  N.  R.  3-STAGE  AXIAL  FLOW  COMPRESSOR 

The  compressor  installation  is  shown  in  Fig.  1.   The  compressor 
inlet  is  in  the  background  of  the  picture.   The  exit  section  is 
equipped  with  a  throttle  valve.   A  torque  meter  is  installed  on  the 
shaft  between  the  drive  motor  and  the  compressor.   The  overall  dimen- 
sions of  the  equipment  are  given  in  Fig.  2. 

A  basic  criterion  of  the  design  of  the  compressor  was  flexibility 
of  operation.   Each  blade  row  may  be  removed.   The  angle  setting  of 
each  blade  is  adjustable  in  0.5  increments.   In  order  to  minimize 
flow  disturbances  by  pressure  probes  the  outside  diameter  of  the  blad- 
ing is  36  in.   For  the  same  reason,,  the  compressor  has  a  hub/tip  ratio 
of  0.6,  giving  an  inside  diameter  of  21.6  in.  and  a  blade  height  of 
7.2  in.   To  permit  simulation  of  conditions  in  a  multistage  unit,  the 
machine  has  three  stages.   One  row  of  inlet  guide  vanes  simulates  the 
effects  of  previous  stages.   Two  rows  of  exit  guide  vanes  remove  the 
whirl  component  after  the  third  stage.   To  insure  that  general  align- 
ment and  tip  clearances  are  maintained  the  compressor  casing  is  very 
rigid. 

The  outer  casing  consists  of  two  cast  iron  half-cylindrical 
shells  bolted  together  along  the  horizontal  plane  through  the  axis. 
The  rows  of  inlet  guide  vanes,  stator  blades,  and  exit  guide  vanes 
are  held  by  bolts  extending  through  the  casing  (Fig.  3).   Inside  the 
casing  a  hollow  steel  shaft  carries  three  cast  iron  drums  for  each 
of  the  rows  of  rotor  blades.   The  shaft  is  supported  by  two  roller 
bearings  whose  outer  races  are  pressed  into  the  supporting  strut 
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assemblies  (Fig.  4).   These  assemblies  rest  in  close  fitting  channels 
in  the  casing  (Fig.  5).   The  six  assembly  struts  have  symmetrical  air- 
foil sections.   Figure  6  shows  the  nine  blade  rows  after  assembly. 

The  casing  was  designed  to  give  maximum  accessibility  for  mea- 
suring instruments.   Six  rectangular  instrument  ports  are  located  in 
the  upper  half  at  30   from  either  side  of  the  vertical  (Fig.  3).   They 
are  placed  in  the  second  and  third  rotor  planes,  the  first,  second  and 
third  stator  planes,  and  behind  the  third  stator.   The  ports  hold  a 
special  instrument  carriage  permitting  detailed  flow  surveys  in  any 
axial  plane.   The  casing  has  numerous  radial  survey  holes  whose  loca- 
tions are  specified  in  Fig.  7. 

The  blades  are  made  of  ALCOA  356  aluminum  alloy  without  heat  treat- 
ment.  There  are  thirty  rotor  blades  and  thirty-two  stator  blades  per 
row.   The  details  of  construction  of  the  blades  are  illustrated  in 
Figs.  8  and  9.   The  blading  is  designed  to  have  a  degree  of  reaction 
of  0.5  at  the  mean  radius.   The  rotor  blades  are  twisted  by  49   from 
hub  to  tip,  and  the  stator  blades  by  13  .   At  the  mean  radius,  the 
important  blade  parameters  are: 

Parameter  Rotor         Stator 

Camber  angle 

Design  stagger  angle 

Thickness-to-chord  ratio 

Chord  length 

The  values  of  the  blade  parameters  at  other  radii,  and  the  method 
used  for  calculating  the  thickness  distribution,  are  specified  in 
Ref.  5. 


20.32° 

29.98° 

43.80° 

28.80° 

0.1 

0.1 

2.60  in. 

2.60  in. 
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The  point  of  maximum  thickness  was  set  at  0.35  of  the  chord 
length  from  the  leading  edge.   The  thickness  was  modified  over  the 
rear  15  per  cent  of  the  section  to  provide  a  trailing  edge  thickness 
0.02  in.   The  thickness  distribution  is  applied  about  a  parabolic 
mean  camber  line.   Figure  10  is  a  picture  of  a  "free-vortex"  stator 
blade.   Tip  clearances  of  0.020  in.  for  stator,  and  0.037  in.  for 
rotor,  are  maintained  by  the  use  of  shims  at  the  point  of  attachment. 

The  attaching  device  permits  variation  of  rotor  and  stator  blade 
stagger  angle.   Figure  11  shows  how  these  changes  are  made.   The  blade 
shaft  extends  through  the  rotor  drum  at  the  hub.   An  adjusting  plate 
is  secured  to  the  blade  shaft  by  a  tapered  pin.   The  plate  has  a  series 
of  eleven  taps  on  a  circle  about  the  shaft  axis  which  are  spaced  at 
intervals  of  4.5  .   Between  the  adjusting  plate  and  the  rotor  drum  is 
a  fixed  sector  which  is  attached  to  the  drum  by  a  pin  which  guarantees 
its  alignment.   The  fixed  sector  has  eleven  holes  which  are  spaced  at 
intervals  of  4.0  of  arc  about  the  blade  shaft  axis.   Blade  alignment 
is  maintained  by  a  set  screw  through  both  holes  in  the  adjusting  plate 
and  the  fixed  sector.   If  the  center  holes  of  the  plate  and  the  sector 
are  lined  up,  the  blade  is  at  the  design  stagger  angle  of  43.8   for 
the  "free-vortex"  rotor.   An  angle  change  of  4.0  can  be  accomplished 
by  keeping  the  set  screw  in  the  center  hole  of  the  adjusting  plate  and 
moving  the  blade  so  the  screw  is  inserted  into  the  next  hole  of  the 
fixed  sector.   Small  angle  changes  of  0.5  are  accomplished  by  moving 
the  blade  slightly  so  that  the  two  holes  immediately  adjacent  to  the 
center  holes  are  lined  up  since  blade  angle  is  changed  only  by  the 
difference  in  arc  between  the  two  holes.   A  similar  arrangement  is  used 
for  the  stator  blade  settings  (Fig.  10).   Figure  12  illustrates  stator 
blade  settings  which  are  one  degree  smaller  than  the  design  stagger 
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angle.   The  two  blades  on  the  right-hand  side  of  the  figure  belong  to 
the  first  stator  row.   The  set  screw  is  two  holes  away  from  the  center 
or  reference  holes.   The  blades  on  the  left-hand  side  of  Fig.  12  belong 
to  the  row  of  inlet  guide  vanes.   They  are  set  at  an  angle  which  is  by 
0.5  smaller  than  the  design  stagger  angle. 

The  assessories  to  the  basic  compressor  will  be  described  by  fol- 
lowing a  path  from  the  inlet  to  the  exit  and  to  the  power  source.   The 
instrumentation  of  the  compressor  will  be  described  in  its  appropriate 
place  in  the  same  sequence.   Details  of  instrument  calibrations  are 
contained  in  Appendix  A. 

The  inlet  duct  is  seen  in  the  background  of  Fig.  1.   It  consists 
of  a  screen,  an  entrance  bellmouth,  and  a  length  of  straight  pipe. 
Figure  13  shows  the  large  mesh  screen  which  prevents  the  ingestion  of 
foreign  matter  from  the  compressor  bay  apron.   Mounted  on  the  screen 
is  a  mercury  thermometer  used  to  determine  ambient  temperature.   It 
has  provisions  for  psychrometric  analysis  of  the  incoming  air« 
Immediately  behind  the  screen  is  the  bellmouth.   It  changes  the  dia- 
meter rather  abruptly  from  almost  two  diameters  at  the  flare,  to 
36  in.  in  the  axial  distance  of  about  15  in.   A  cylindrical  inlet  duct 
two  diameters  long  connects  the  bellmouth  to  the  compressor  proper. 
The  inner  surface  of  the  duct  is  enameled  to  give  smooth  flow  surfaces. 

The  flow  rate  through  the  compressor  was  determined  by  surveys  in 
the  inlet  duct.   Details  of  the  calibration  are  given  in  Appendix  B. 
The  survey  plane  is  about  midway  between  the  bellmouth  and  the  com- 
pressor.  The  surveys  were  made  with  a  Prandtl-type  Pitot  tube.   In 
Fig.  14  the  probe  is  installed  for  taking  a  vertical  traverse.   It  is 
possible  to  disassemble  the  probe  while  the  compressor  is  running  (Fig.  15). 
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It  may  then  be  remounted  for  a  horizontal  traverse.   Traverses  at  30 
and  60   from  the  vertical  direction  are  possible.   These  traverses  can 
be  made  by  unbolting  the  entire  inlet  duct  from  the  compressor  and 
rotating  it  on  its  cradle.   The  design  of  the  probe  prevents  measure- 
ments closer  than  0.25  in.  to  the  inlet  duct  wall. 

A  honeycomb  straightener  is  installed  behind  the  inlet  Pitot  tube 
to  equalize  the  flow  entering  the  compressor  (Fig.  14).   Attached  to 
the  forward  strut  assembly  is  an  ogival  wooden  fairing  (Fig.  5)  to 
provide  a  smooth  transition  of  the  flow  from  the  inlet  duct  to  the 
compressor  annulus. 

Pressure  measurements  may  be  taken  at  any  one  of  the  radial 
survey  ports  described.   Figure  16  shows  the  probe  used  for  these 
measurements.   It  is  a  United  Sensor  three-hole  probe,  model  YC-120. 
A  central  hole  measures  stagnation  or  total  pressure.   A  static  port  is 
located  on  each  of  two  faces  of  a  wedge  at  an  angle  of  45  from  the 
total  pressure  hole.   The  static  pressures  are  balanced  on  a  water 
manometer  to  insure  alignment  of  the  probe  in  flow  direction.   The 
probe  is  mounted  on  a  probe  holder  which  has  vernier  scales  permitting 
radius  adjustments  to  0.01  in.  and  angle  adjustments  to  0.1  . 

The  accurate  determination  of  flow  rate  is  essential  to  compressor 
analysis.   For  this  purpose  a  modified  Prandtl  type  Pitot-static  probe 
was  installed.   Figure  17  shows  the  installation  details.   The  probe 
is  inserted  in  a  radial  survey  hole  located  ahead  of  the  inlet  guide 
vane  row.   It  protrudes  forward  into  a  channel  between  two  adjacent 
struts  of  the  supporting  assembly  which  is  shown  in  Fig.  5.   The  probe 
is  approximately  at  the  center  of  the  channel  between  struts.   The 
probe  is  positioned  in  the  radial  survey  hole  by  a  brass  plug,  and 
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secured  to  the  casing  by  a  simple  locking  device.   Figure  18  shows  the 
probe  and  its  locking  device  in  the  background.   In  the  middle  of  the 
figure  are  seen  the  three-hole  probe  and  its  holder  located  between  the 
inlet  guide  vanes  and  the  first  rotor.   In  the  foreground,  a  plug, 
originally  located  at  the  same  axial  position  as  the  probe,  has  been 
removed  to  show  a  hole  through  which  the  probe  can  be  inserted. 

Pressures  obtained  by  the  probes  were  measured  by  a  Texas  Instru- 
ments Fused  Quartz  Precision  Pressure  Gauge  which  is  shown  in  Fig.  19. 
The  velocity  head  in  the  inlet  is  about  0»6  in.  of  water.   A  low- 
pressure  bourdon  tube  was  obtained  for  the  gauge  which  measures  from 
0  to  166  in.  of  water  with  the  instrument  reading  from  0  to  200,000 
counts.   The  calibration  constant  of  the  low  pressure  bourdon  tube  was 
determined  to  be  240*423  counts/psf  of  the  pressure  gauge  (See  Appendix 
A.)«   The  pressure  gauge  operates  both  in  manual  and  servo  modes,  the 
latter  providing  automatic  nulling  of  the  unit. 

Since  flexibility  of  pressure  selection  was  desirable,  two 
Giannini  Sp-lOlA  pressure  scanners  were  connected  to  the  pressure  and 
reference  ports  of  the  bourdon  tube.   One  of  the  scanners  is  seen  in 
the  foreground  of  Fig.  19.   The  arrangement  of  pressures  to  the  12- 
channel  switches  is  shown  in  Fig.  20.   Later  in  the  study  a  simpler 
system  of  five  valves  and  two  manifolds  was  constructed  (Fig.  21)  to 
eliminate  leakage  flows  that  seem  to  have  occurred  in  the  scanners. 
At  this  time  the  original  3-hole  probe  was  placed  permanently  behind 
the  third  stator  (Fig.  22).   Another  3-hole  probe  was  placed  ahead  of 
the  first  stage  rotor.   This  arrangement  permitted  direct  measurement 
of  the  pressure  increase  across  the  three  stages. A  P»   . 
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Figure  22  also  shows  the  instrument  traverse  carriage  mentioned 
before.   It  was  used  for  a  survey  of  total  pressure  in  peripheral 
direction  behind  the  third  stator  at  two  different  radii  (See  Appendix 
C.)«   It  i-s  possible  to  vary  the  radial  locations  of  the  probe  with  an 

accuracy  of  0.01  in.   The  probe  may  be  rotated  through  360  about  its 

o 
axis.   The  carriage  can  move  the  probe  15  in  peripheral  direction, 

which  covers  a  whole  blade  spacing  since  a  stator  blade  channel  covers 

11.25  ,  and  a  rotor  blade  channel  12.0   in  peripheral  direction.   The 

locations  of  the  radial  survey  holes  with  reference  to  the  indicated 

meridional  angle  of  the  carriage  are  shown  in  Fig.  23.   The  represented 

stator  profiles  are  at  the  design  stagger  angle  at  the  mean  radius. 

The  actual  stagger  angle  measured  with  respect  to  a  plane  through  the 

axis  is  the  complement  of  the  indicated  angle. 

A  short  cylindrical  duct  connects  the  compressor  to  the  exit 
elbow, which  is  equipped  with  sheet  metal  turning  vanes  to  minimize 
losses.   A  transition  piece  connects  the  elbow  to  the  throttle  valve 
which  is  seen  in  Fig.  1.   The  valve  consists  of  two  rectangular  metal 
doors  which  move  in  slides.   The  position  of  the  doors  is  controlled  by 
a  right-  and  left-hand  threaded  lead  screw  which  is  rotated  by  a 
Graham  variable-speed  transmission.   A  revolution  counter  indicates 
the  approximate  valve  opening. 

The  compressor  is  driven  by  a  50  HP  Fairbanks  Morse  induction 
motor.   It  requires  a  three-phase,  sixty-cycle,  440  volt  power  supply, 
and  can  operate  at  two  fixed  speeds;  namely,  at  about  900  rpm  or 
1200  rpm. 

Attached  to  the  motor  shaft  is  a  Baldwin-Lima-Hamilton  SR-4 
torque  meter  type  A.   It  is  shown  in  Fig.  24  during  calibration  with 
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known  weights  attached  to  a  lever.   In  the  background  may  be  seen  a 
Brown  Instruments  strain  gauge  readout  which  is  of  the  standard 
Wheatstone  bridge  type.   A  conversion  constant  of  4.1565  ft-lb/lb  of 
torque  meter  readout  was  determined  (See  Appendix  A.).   Under  the 
protective  cover  next  to  the  torque  meter  in  Fig.  24  is  installed  a 
six  lobe  flux  cutter  to  measure  the  motor:  speed  by  means  of  an 
electronic  counter. 
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SECTION  3 
FLOW  RATE  CALIBRATION 

The  permanent  Pitot-static  tube  was  installed  to  provide  rapid 
and  accurate  measurements  of  volume  flow  rate.   At  a  given  throttle 
valve  setting  this  probe  was  calibrated  against  the  result  of  two 
traverses  of  the  cylindrical  inlet  duct.   Through-flow  velocities  in 
the  flow  annulus  ahead  of  the  inlet  guide  vanes  were  varied  from 
110  ft/sec  to  70  ft/sec,  so  that  compressibility  effects  could  be 
ignored.   The  objective  of  the  calibration  was  to  determine  a  cali- 
bration constant  which  gave  the  volume  flow  rate  when  multiplied  by 
the  velocity  measured  at  the  permanent  probe. 

An  initial  survey  of  the  inlet  duct  was  made  at  the  highest  pos- 

3 
sible  flow  rate  of  512  ft  /sec.   Figure  25  shows  a  plot  of  the  mea- 
sured velocities  of  both  horizontal  and  vertical  traverses.   The 
smooth  velocity  profile  in  the  duct  is  evident.   Figure  26  is  a  graph 

of  the  measured  velocities  in  the  inlet  duct  at  the  lower  flow  rate  of 

3 
310  ft  /sec,  which  shows  a  higher  degree  of  scatter.   For  these  tests 

the  pressure  indicator  was  operated  in  the  servo  mode.   Even  in  the 

meter  mode  the  indicated  pressure  oscillated  considerably.   Pressure 

fluctuations  seemed  to  have  been  initiated  by  the  vibration  of  the 

brass  rod  on  which  the  traversing  Pitot-static  tube  was  mounted  in  the 

inlet  duct.   The  pressure  indicator  corrects  a  disparity  between  the 

angular  position  of  a  mirror  of  the  bourdon  tube  and  the  indicator 

dial  at  the  rate  of  the  full  scale  reading  of  200,000  counts  in  120 

seconds.   With  the  oscillating  pressure  applied  to  the  sensitive 

instrument,  oscillations  of  50  counts  were  observed  for  a  pressure 
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corresponding  to  a  reading  of  600  counts.   To  reduce  these  oscillations, 
variable  length  capillary  tube  damping  devices  were  installed  both  in 
the  total  and  static  pressure  lines  of  the  traverse  Pitot-static  tube. 
Details  of  the  arrangement  and  the  results  of  tests  are  shown  in 
Appendix  B-l. 

A  special  survey  was  made  to  determine  the  thickness  of  the  bound- 
ary layer  on  the  walls  of  the  cylindrical  inlet  duct.   However,  con- 
struction of  the  traverse  probe  prevented  taking  readings  closer  than 
0.25  in.  from  the  wall.   Figure  27  is  a  plot  of  the  readings  at  high 
and  low  flow  rates.   In  each  case  it  is  evident  that  large  losses  occur 
near  the  walls.   With  the  velocity  distributions  of  Fig.  27  the  flow 
rate  is  about  2.5  per  cent  smaller  than  the  value  obtained  without 
considering  the  changes  near  the  walls . 

Data  were  taken  in  the  boundary  layer  at  intervals  of  0.25  in. 
from  17.75  in.  to  16.00  in.   In  the  main  stream, data  were  taken  at 
4.0  in.  intervals  from  16.00  in.  to  the  centerline.   With  horizontal 
and  vertical  traverses,  four  points  were  obtained  at  each  radius  and 
two  at  the  axis  for  a  total  of  46  locations.   At  each  station  the 
following  measurements  were  taken: 

q.     -  velocity  head  in  the  inlet  duct  (counts) 

q     -  velocity  head  of  the  permanent  Pitot-static  probe  (counts) 

P  .    -  static  pressure  with  reference  to  atmosphere  in  inlet 
si  r 

duct  (counts) 
P     -  static  pressure  with  reference  to  atmosphere  at  permanent 

Pitot-static   probe    (counts) 

T     -  temperature  at  the  inlet  screen  (  F) 

P     -  barometric  pressure  (in.  of  Hg) 
Ba  r 

T     -  temperature  of  column  of  mercury  (  F) 
Ba  r 
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The  barometric  pressure  P    (in.  Hg)  was  converted  into  pounds 

Bar 

per  square  foot  to  obtain  the  atmospheric  pressure  P  by 

A 

PA  -   P8ar   (71.  4G7)  (Psfa)  (1) 

where  the  constant  was  obtained  by  correcting  the  height  of  the  column 

of  mercury  for  temperature  variations  in  specific  gravity  of  the  liquid. 

The  inlet  temperature  was  converted  to  absolute  temperature, and  was 

considered  to  be  the  total  temperature  T   because  of  the  low  veloc- 

r         to 

ities,  or 

T,D  =  \     -  4^7  (°R)  (2) 

The  static  pressure  P  at  the  inlet  was  converted  from  the  meter 
reading  to  an  absolute  pressure  by 

A     £4-0.42.3 
where  the  constant  was  obtained  by  the  calibration  procedure  described 

in  Appendix  A.   The  pressure  P  measured  at  the  permanent  probe  was 

obtained  in  a  similar  manner.   The  local  air  density  p.  at  the  inlet 

was  computed  from 


_  O.OOZ^aPc   W  =  lL^xl0-)  (lb  sec2/ft4)      (4) 

and  similarly  p   at  the  permanent  probe.   The  velocity  head  in  the 
inlet  q..  was  converted  from  the  meter  reading, 

q    =   ^ (psf)  (5) 

b>1    240.  413 

The  dynamic  head  at  the  permanent  probe  q   was  obtained  similarly. 

The  velocity  V.  in  the  inlet  is  then 
J      i 


Vt  -V  ^ o%°X  (ft/sec)  (6) 
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A  similar  procedure  gave  the  velocity  V  at  the  permanent  probe.   The 

arithmetic  average  of  the  46  values  of  V  obtained  during  the  carrying 

out  of  the  traverses  gives  the  average  value  V    of 

°  oav 

Voa^  =  J£     Z-i  V°j  (ft/sec) 

The  measured  values  of  V.  were  adjusted  to  account  for  variations  in 
flow  rate  due  to  slight  compressor  speed  variations 

Visd  "  Vt  — (ft/sec) 

The  four  values  of  V.  ,  obtained  at  each  radius  were  averaged  by 

4- 

V.av  =  "5"  Z-  V^j  (ft/sec) 

J-1 
The  volume  flow  rate  was  first  calculated  by  neglecting  the  losses  in 

the  boundary  layer  to  give  N^ j.  .   In  general  the  volume  flow  rate 

is  obtained  by 

V  r  Jr  (ft3/sec) 

If  the  above  equation  were  integrated  directly,  the  velocity  at  the 
axis  would  not  be  included  in  the  volume  flow  rate  calculation 
because  the  zero  radius  is  a  member  of  the  product  under  the  integral. 
Therefore  the  integral  was  replaced  by  a  summation. 

\Ft   =  H     V.'av   *A  (ft3/sec) 

The  quantity  A  A  =  2  Tf  r  h/144  represented  that  part  of  the  annulus 

area  in  the  center  of  which  V,    was  measured,  where  h  is  the  height 

iav 

of  the  area.   This  method  assumed  that  the  velocity  measured  at  the 
radius  of  16.0  in.  exists  also  at  the  wall.   A  volume  flow  rate 
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deficiency  was  then  calculated  to  account  for  the  difference  in 
velocities  at  16.0  in.  radius  and  those  measured  in  the  boundary 
layer,  or 

VD  *   Vtav  ,fc.  -  VCaV  9>L-  (ft/sec) 

By  a  trapezoidal  integration  the  volume  flow  deficiency  vV  becomes 

,  Jfc.O 

\F0    s      I  7f      I  VD       r     Jr  (ft3/sec) 

7  /fc.O 

Then  the  actual  volume  flow  rate  is 

^"c  -   ^"t    ^p  (ft3/sec) 

A  calibration  constant  CCC  for  each  run  was  obtained  from 

^C  3 

CCC  =  — - — — (ft  /sec  per  ft/sec) 

oav 
A  first  series  of  calibration  runs  was  invalidated  when  a  leakage 

hole  was  discovered  in  the  traversing  Pitot-static  tube.   After  repairs, 

the  calibration  constants  CCC  were  obtained  with  a  maximum  relative 

error  of  0.1  per  cent.   The  final  calibration  constant  was  4.4050 

3 
ft  /sec  per  ft/sec  velocity  obtained  from  the  readings  of  the  perma- 
nently installed  Pitot-static  probe. 

A  data  reduction  computer  program,  ONRFLO,  for  these  calibration 
procedures^ was  generated  for  the  I.  B.  M.  Model  360  computer.   A 
listing  of  the  program  is  included  as  Appendix  B-2.   Sample  output 
data  from  the  final  run  is  given  in  Appendix  B-3. 
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SECTION  4 
MEASUREMENT  OF  COMPRESSOR  PERFORMANCE 

The  objective  of  the  study  was  to  determine  the  effect  of  changes 
of  the  stator  blade  orientation  on  compressor  performance.   It  was 
decided  to  restrict  performance  investigation  to  the  three  stages  only; 
namely,  from  a  station  ahead  of  the  first  rotor  to  a  station  after  the 
last  stator,  or  between  locations  SP-2  and  SP-8  of  Fig..  7*   The  effects 
of  the  inlet,  the  inlet  guide  vanes,  the  exit  guide  vanes,  and  the 
diffuser  were  not  considered. 

For  tests  at  a  particular  stator  stagger  angle  the  barometric 

pressure  P    was  determined  first.   For  various  settings  of  the 
Bar 

throttle  valve  the  following  measurements  were  taken: 

T     -  ambient  temperature  (  F) 
o 

P     -  static  pressure  with  reference  to  atmosphere  at  the 

permanent  Pitot-static  probe  (counts) 
q     -  velocity  head  at  the  permanent  Pitot-static  probe  (counts) 
P  .    -  total  pressure  with  reference  to  atmospheric  between  the 

inlet  guide  vanes  and  the  first  rotor  (counts) 
A  P    -  total  pressure  difference  across  the  three  stages  (counts) 
T     -  reading  of  torque  meter  (lb) 

N     -  compressor  speed  measured  by  electronic  counter  (rpm) 
With  the  relations  listed  on  p.  35  these  measuring  data  establish 
the  following  quantities: 

P*   =  ?b&r     C-7/.4C7;  (psfa)  (1) 

Tt.  '  T0  *  4  5^.7  (°r)  (2) 

P   -    p^° 

Ko  240   4 11  <psfa>  <3> 
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=  0.00^78   PQ   5,8.7  .    _?^  ^U^)  (lb  Sec2/ft4)  (4) 

(»Utt)(l44)      T\0  Ti0 


bo1         240.4£3 


V.     -J   ^^ 


(ft/sec)  (6) 


The  absolute  total  pressure  ahead  of  the  first  stage,  or  after  the 

inlet  guide  vane,  P  .   ,  is 

tiga' 


p.     .   p.  -  Ln 


"9«         A       i^Til  (ps£)  (7) 

The  pressure  difference  across  the  three  stages  A  P     is  obtained 
from 


aP 


»•«*  =  210. -Ui  <psfa) 


From  the  torque  meter  reading  T   the  actual  torque  T   is 

T    -    Tr  (4. 1 5-  ts>         (ft-ibf) 

where  the  constant  was  obtained  by  the  calibration  procedure  described 
in  Appendix  A.   The  angular  velocity  CO  is  with  the  measured  speed,  N 

ZTf  N 
CO  b   — (rad/sec) 

4,0 

Both  volume  and  weight  flow  rates  were  calculated.   The  volume 
flow  rate   v  ^   is  then 

\fc      -   V0   (4.4050)  (ft3/sec)  (8) 

where  the  constant  was  obtained  by  the  calibration  procedure  described 
in  the  preceding  section  and  Appendix  B.   The  weight  flow  rate  w  is 


39 


.      .p    32.H4   P*  518.7       .c-    P»    ;.,  „„,, 

w  =  c  o*wxi«)Tlo  =  ^  17, U 88fe)    <lbm/sec)  <9) 


For  comparison  with  the  prediction  program  a  referred  flow  rate  w  is 


defined  by 


Wr    =      o • (lbm"V°R/sec   psia)    (10) 

PA  /»44 

Two  basic  indices  of  performance  were  computed;  namely,  the 

3-stage  pressure  ratio  P  -   where 
°   r  r3-s> 

rr  3„s  -  — (ii) 

and   the   3-stage   total-to-total  efficiency   n  ,_<.    obtained   from 

P 
rWfh 

13-5        =  r  (12) 

r  WA 

where  P..,,  and  P„  are  the  theoretical  and  actual  horsepower,  respec- 
Wth      Wa  r  >  r 

tively.   The  theoretical  power  is  P  ,  ,given  by 

PWfh  =  w   A^isn78.3)  (ft-lbf/sec) 

where  Ah.   is  the  isentropic  specific  enthalpy  difference  for  a 

frictionless  compression  from  P^,    to  P  ,  after  the  third  stator, 
r  tiga     td 

Assuming  a  perfect  gas  there  is 


Ak;s  -   cp  A   1*3.3  (Btu/lbm,  °R) 

where  AT  '    is  the  isentropic  temperature  difference  from  the  total 

temperature  T  .    ahead  of  the  first  stage,  to  T  .  after  the  third 
tiga  °       td 

stator  (See  Fig.  28.)-   Since  T  .    very  nearly  equals  the  ambient 

tiga 

temperature  T      ,    there   is 
to 
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K  =  cp  Tto 


/  r  ' 


ta 


With  the  isentropic  relation 


1 


(Btu/lbm,   R) 


r, 


td 


-td 


T 


p. 


the  final  expression  for  the  determination  of  the  theoretical  power  is 


=  W  C. 


'W+h    vv  "\> 
The  actual  power  is  obtained  from 


\0  [(  Pr3.s)  y-  IJ77S.3       (ft-lbf/sec) 


(13) 


»WA     =      Ipq     U)  (ft-lbf/sec) 

The   total- to- total  efficiency  at   the   three  stages   is   then 


a 


3-S 


\v  cp  lt< 


[ttvS)?-  l] 


77  8.3 


(14) 


For  comparison  with  other  compressor  data  certain  non-dimensional 

parameters  are  introduced.   The  average  flow  function  Toy  *s  fc^e 

ratio  of  through  flow  velocity  V  and  peripheral  speed  at  the  mean 

diameter  U   ,  or  with 
av 

V*    =      451.4/144       =    "e   l*«0         (n/sec) 
and,    for   the  mean  diameter  of  14.4   in. 


U. 


ftV 


CO 


14,4 


(ft/sec) 


there   is 


<t> 


AV 


V^c      (i2)(K4) 

CO      (l4.4-)C4S|.4) 


NTc 


(JO 


(0.184Z)  (15) 
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An  average  dimensionless  stage  pressure  coefficient  7  Av/  is  defined  by 


<XV 


^^  =  vl  e6  Uw z       &.f) 


01V 


or 


_A_P; 

3  e.  Uw 


T*v        ~       ~  TT2"  (16) 


The  so-called  dimensionless  deflection  coefficient  L  As/  is  defined  by 


aP  a  P         W  t  t  * 

Al  3-5+i,  ..  &  r^A  _  _^v  p  (J      (psf) 

3       ^?J-S    W 


=  ^  Po  Uj  (p-o 


Hence, 


13-S 

The  performance  calculations  were  made  for  the  conditions  along 
the  mean  radius  of  the  stage, 14. 4  in.,  by  assuming  that  these  conditions 
are  representative  of  the  performance  of  the  3-stage  compressor. 
Detailed  investigations  which  support  this  simplification  are  given  in 
Appendix  C.   The  data  reduction  was  accomplished  by  the  computer 
program  ONRETAs  which  is  listed  also  in  Appendix  C. 

Data  were  taken  on  11  occasions  for  the  following  settings  of 
stator  stagger  angle: 
Run 


1,  2, 

3 

4,  5 

6,  7 

8,  9 

10 

11 

Stator 

Stagger 

Angle 

27, 

,8° 

23. 

.8° 

31. 

,8° 

35, 

,8° 

39, 

,8° 
„o 

44.3 
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SECTION  5 
PERFORMANCE  PREDICTION  PROGRAM 

The  program  used  to  predict  the  performance  of  multi-stage  com- 
pressors is  called  AXC03.   It  is  listed  in  Appendix  D  and  represents 
an  adaptation  of  a  program  by  Vavra  (3),  modified  by  Gibbons  and 
Bartels  (4). 

The  method  assumes  that  the  flow  is  axisymmetric.   Computations 
are  made  for  the  streamlines  at  the  mean  radius  of  a  stage.  A  stage 
by  stage  analysis  of  the  machine  is  performed.   The  output  of  one 
stage  is  used  as  input  for  the  next  stage.   The  computation  is  started 
at  a  flow  rate  less  than  that  anticipated  for  surge,  and  the  flow  rate 
is  then  increased  by  small  increments  to  cover  the  whole  operating 
range  of  the  machine „ 

Program  AXC03  uses  two  principal  subroutines;  namely,  one  called 
ROTOR  and  the  other  called  STATOR.   They  establish  the  velocity 
triangles  of  a  rotor  or  a  stator  row  of  blades,  respectively,  from  the 
known  geometry  of  the  bladings,  flow  rate,  and  discharge  conditions  of 
the  preceding  row  of  blades.   Both  subroutines  make  use  of  three 
additional  subroutines. 

Subroutine  THEORY  calcualates  the  flow  conditions  for  minimum 
profile  losses.   Cascade  data  by  Lieblein  are  used  for  this  purpose  (2). 
These  data  are  presented  as  curves  for  zero-camber  incidence  angle 
(kQ)     ,  slope  factor  n,  zero  camber  deviation  angle  (€0)  .   ,  and 
slope  factor  m,  all  given  as  functions  of  the  inlet  air  angles   6 
or  ai   .  ,  with  the  blading  solidity  J*  as  a  parameter.   The  solidity 
QT    is  defined  as  the  ratio  of  blade  chord  c  and  blade  spacing  s. 
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The  curves  described  above  have  been  expressed  by  polynomials  in  powers 
of  solidity  and  inlet  air  angle  which  are  used  in  subroutine  THEORY. 
The  accuracy  of  the  polynomials  has  been  verified,  and  the  errors  are 
listed  in  Appendix  D.   The  incidence  angle  for  minimum  loss  is  deter- 
mined by  an  iteration  procedure  which  changes  the  inlet  air  angle  until 
agreement  is  reached  within  a  tolerance  of  0.02  .   The  deflection  which 
occurs  at  the  minimum  profile  loss  can  then  be  calculated  with  the 
established  polynomials  also. 

The  subordinate  subroutine  CASCAD  calculates  the  performance  of 
the  blading  at  the  actual  incidence  or  inlet  flow  angle  of  the  blade 
row  by  using  a  graph  of  Lieblein  that  establishes  the  change  of  the 
deviation  angle  with  incidence  angle  ( a  b    fd  C )  as  a  function  of 
inlet  flow  angle  and  solidity. 

In  the  earlier  computer  programs  the  curves  of  this  graph  were 
expressed  by  a  polynomial  that  covered  the  range  of  solidities  from 
0.0  to  1.8.  Since  the  curves  have  nearly  exponential  character,  which 
cannot  be  expressed  with  ease  by  polynomial  functions,  it  was  found 
that  an  error  of  the  order  of  10  per  cent  was  possible.  Therefore, 
a  new  polynomial  was  established  for  solidities  between  0.6  and  1.4 
which  limits  the  error  to  less  than  1.9  per  cent. 

Subroutine  CASCAD  further  utilizes  experimental  limits  on  blading 

performance  by  Howell  (9)  which  are  illustrated  in  Fig.  29.   Howell 

found  that  for  most  cascades  the  ratio  of  the  profile  drag  coefficient 

C_.  at  an  arbitrary  incidence  angle  L     and  C_  .   at  the  minimum  loss 
D  °  Cmin 

incidence  L t0    is  a  unique  function  of  the  quantity  (L"t,0)/6|0  , 

where  6  |0  is  the  optimum  flow  deflection  for  the  incidence  angle  I  ( Q . 

Since  the  applicability  of  the  data  presented  in  Ref.  2  is  limited  to 
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design  point  calculations,  where  the  incidence  angles  are  close  to  the 
nominal  incidence  Lq,  the  curves  cannot  be  used  for  off-design  analy- 
ses without  restrictions.   Without  such  limits  the  flow  deflections  € 
could  be  increased  simply  by  increasing  the  incidence  angle  L  y  since 
the  data  of  Ref.  2  do  not  establish  criteria  for  flow  separations 
which  are  associated  with  a  radical  increase  in  profile  losses.  Exper- 
imental results  of  Howell  show  that  the  flow  deflection  6    at  incidence 
C  should  not  exceed  the  optimum  deflection  £  l0  by  more  than  25  per  cent 
to  avoid  excessive  losses.   He  states  that  this  limit  occurs  if  the 
profile  loss  coefficient  C_  is  about  twice  the  minimum  loss  coefficient 
C_  .  .   The  solid  curve  in  Fig.  29,  labeled  e/e|0,  is  supposed  to  be 
the  relation  between  £/t/0  an^  (£" L  t0)/ £  ,0    as  obtained  from  subroutine 
CASCAD  with  the  data  from  Ref.  2.   It  will  be  assumed  that  this  curve 
holds  for  negative  values  of  (i"l|*)/6..   only.   For  positive  values 
of  this  parameter  the  dashed  curve  will  be  applied  to  satisfy  the 
limitations  of  Howell.   This  curve  is  somewhat  arbitrary  but  represents 
the  best  estimate  for  the  actual  operating  condition  of  compressor 

cascades.   Figure  29  further  shows  a  curve  relating  the  ratio  C^/C^  . 

•  D  Dmm 

to  (L~ L j0) I  €  lo     ,  which  is  adapted  from  Ref.  9.   Both  this  curve  and 
the  dashed  curve  for  the  deflection  ratio  have  been  expressed  analy- 
tically in  subroutine  CASCAD.   If  the  ratio  C^/C^  .   exceeds  a  value  of 
J  D  Dmin 

2.0  for  positive  quantities  (l-l(0  )/fe(0  >  the  program  prints  out 

"surge  in  blade  row  xx"  and  proceeds  with  calculations  at  an  increased 

flow  rate.   The  positive  value  of  ( C~  ^,0)/£  ._  at  which  this  occurs  is 

the  input  quantity  SMAX.   If  C^/C^_.   exceeds  a  value  of  2.0  at 
r        n      J  D  Dmin 

incidence  angles  L    smaller  than  C  (  0  ,  the  program  prints  out  "minimum 
pressure  in  blade  row  xx"  and  halts  computation.   The  negative  value  of 
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(C-L|0)/&10   at  which  this  occurs  is  the  input  quantity  SMI.   In 

both  instances  the  symbol  xx  refers  to  the  number  of  the  blade  row, 

and  indications  are  given  also  whether  the  row  is  a  stator  or  a  rotor. 

For  any  incidence  angle  within  the  useful  operating  range, the  values 

of  the  deflection  ratio  and  of  the  ratio  C^/C^  .   are  computed  by  sub- 

D  Dmin  J 

routine  CASCAD  and  control  returned  to  the  blade  row  subroutines 
ROTOR  or  STATOR. 

The  efficiency  of  a  blade  row  is  determined  by  subroutine  ETACAL. 
Losses  due  to  tip  clearances,  secondary  flows,  and  wall  friction  are 
taken  into  account  with  the  relations  proposed  by  Vavra  (1), 

If  conditions  of  surge  exist,  the  flow  rate  is  increased  by  about 
1.0  per  cent  and  the  entire  computation  is  repeated.   If  the  incidence 
angles  are  in  the  so-called  useful  operating  region  of  Fig.  29,  the 
output  of  subroutine  ROTOR  or  STATOR  becomes  the  input  for  the  sub- 
routine that  calculates  the  next  row  of  blades.   Interstage  data  are 
printed  out  as  the  program  progresses  from  one  stage  to  the  next.   If 
all  blade  rows  of  a  machine  have  been  processed,  the  overall  compressor 
efficiency  and  pressure  ratio  are  computed  for  the  particular  flow  rate. 
This  flow  rate  is  then  increased  by  a  specified  amount,  and  the  calcu- 
lating process  is  repeated.   If  the  so-called  minimum  pressure  is 
reached,  the  computation  stops  and  a  summary  of  the  overall  performance 
parameters  is  printed  out. 

The  details  of  the  calculation  of  the  flow  through  the  inlet  duct, 
the  inlet  guide  vanes,  the  exit  guide  vanes,  the  diffuser  and  the  dis- 
charge passages  are  not  described  in  this  thesis  since  it  is  concerned 
only  with  the  performance  prediction  of  the  three  stages  of  the  com- 
pressor from  a  station  ahead  of  the  first  rotor  to  the  discharge  at 
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the  third  stator.   The  methods  applied  for  the  flow  analysis  in  these 
passages  may  be  found  in  Vavra  (3)  and  Gibbons  and  Bartels  (4). 

According  to  Howell  there  occur  additional  effects  that  influence 
the  off-design  calculations.   The  growth  of  boundary  layers  on  the 
walls  of  the  annular  flow  channel  will  change  the  velocity  profiles 
and  can  reduce  significantly  the  effective  flow  area.   The  program 
takes  account  of  this  effect  by  a  blockage  factor  that  can  vary  from 
stage  to  stage.   This  boundary  layer  growth  is  responsible  also  for 
increasing  peaks  in  the  radial  distribution  of  the  axial  velocity 
components  in  succesive  stages.   These  higher  velocities  outside  of 
the  wall  boundary  layers  produce  smaller  incidence  angles  for  the 
main  portion  of  the  flow,  thereby  decreasing  the  actual  flow  deflection 
and  reducing  the  work  absorbed  by  the  fluid  in  the  stage.   Hence 
in  a  compressor  that  consists  of  stages  with  identical  bladings, 
the  last  stages  will  produce  smaller  pressure  ratios  than  the  stages 
at  the  compressor  inlet.   Experience  shows  the  efficiency  in  successive 
stages  is  not  reduced  by  the  peaking  of  the  velocity  profiles,  and 
their  effects  are  usually  taken  into  account  by  a  so-called  work-done 
factor  that  is  about  unity  for  the  first  stage  and  gradually  decreases 
for  the  successive  stages.   Such  work-done  factors  can  be  introduced 
in  program  AXC03,  but  for  the  present  compressor  they  will  be  taken  as 
unity  because  of  the  small  number  of  stages  and  the  large  blade 
heights.   For  the  same  reasons  the  blockage  factor  is  assumed  to  be 
equal  to  unity  also. 

The  geometry  of  the  blades  could  be  obtained  from  the  information 
of  Ref.  5.   Values  at  the  mean  radius  were  used  for  the  analysis. 
Bowen  (5)  determined  that  the  average  flow  angle  at  the  mean  radius 
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after  the  inlet  guide  vanes  ALE  was  20.0   for  the  design  inlet  guide 
vane  stagger  angle  of  11.4  ,  independent  of  flow  rate.   From  Fig.  12 

it  can  be  recognized  that  the  inlet  guide  vanes  were  set  at  a  stagger 

o  o 

angle  of  10.9   for  the  tests.   Therefore  a  value  of  19.5  was  used  in 

the  analysis  program  for  the  average  flow  angle  leaving  the  inlet 

guide  vanes.   Reference  5  does  not  give  values  of  C__.   for  the 

Dm  in 

blading.   It  is  possible  to  introduce  in  the  program  various  values  of 

C__.   both  for  rotor  and  stator  blading.   Hence  it  was  possible  to 
Dmin  °  r 

estimate  an  average  apparent  value  of  C_  .   for  blading  by  comparing 

Dm  in 

measured  results  to  predictions  of  the  analysis  program  for  several 

assumed  values  of  Cn  .  . 

Dm  in 

A  prediction  of  the  3-stage  efficiency  for  different  stator 

stagger  angles  is  shown  in  Fig.  30  for  an  assumed  value  of  C_  .   of 

Dm  in 

0.020.   The  light  line  connecting  the  individual  peak  efficiencies 

represents  the  calculated  operating  envelope  which  could  be  obtained 

if  continuous  control  of  stator  stagger  angle  for  maximum  efficiency 

were  possible. 

For  comparison  with  measured  performance  a  series  of  calculations 

were  made  for  stator  stagger  angles  of  23.8  ,  27.8  ,  31.8  ,  35.8  , 

39.8  ,  and  44-3  .   For  each  stagger  angle,  values  of  C^  .   of  0.000, 

Dm  in 

0.006,  and  0.008  were  assumed.   To  avoid  confusion  with  the  computa- 
tions of  actual  performance  based  on  measured  quantities,  which  have 
been  called  runs,  each  calculation  of  the  analytical  prediction 
program  has  been  assigned  a  two-digit  case  number.   The  first  digit 
identifies  the  stator  stagger  angle  of  the  case  as  follows: 
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First  Digit  Stator  Stagger  Angle 

1  23.8° 

2  27.8° 

3  31.8° 

4  35.8° 

5  39 . 8° 

6  44.3° 

The  second  digit  identifies  the  assumed  value  of  C^  .   of  the  case 

Dmin 

as  follows: 

,  ^.  .                Assumed  Value  of  C^  . 
Second  Digit  Dmin 

1  0.000 

2  0.006 

3  0.008 

For  example,  case  42  is  the  predicted  performance  at  a  stator  stagger 

angle  of  35.8  with  a  value  of  (L  .   of  0.006.   The  summary  output  of 

Dmin 


each  case  is  presented  in  Appendix  D. 
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SECTION  6 
DISCUSSION  OF  RESULTS 

Experimental  data  for  performance  measurements  were  taken 
during  55  hours  of  running  time.   Another  81  hours  of  operation  were 
used  for  miscellaneous  calibration  and  testing,  including  44  hours 
for  flow  rate  calibrations. 

The  rotor  stagger  angle  was  set  at  43.8  for  all  runs.   The 
design  stator  stagger  angle  is  288  .   For  the  performance  tests,  the 
stator  stagger  angle  was  set  at  23.8°,  27.8°,  31.8°,  35.8°,  39.8°, 
and  44.3  .   Measured  performance  is  presented  in  comparison  with  the 
predicted  results  of  AXC03,and  also  by  establishing  dimensionless 
performance  parameters.   Figures  31  through  36  are  graphs  of  measured 
performance  in  comparison  with  the  prediction  by  AXC03  for  each  stator 
blade  angle  setting.   Figures  37  and  38  are  summary  plots  of  maximum 
efficiency  and  pressure  ratio  for  all  stator  blade  angle  settings. 
Figures  39  through  44  are  graphs  of  the  non-dimensional  parameters  at 
each  stator  stagger  angle.   Figure  45  is  a  summary  plot  of  maximum 
efficiencies  and  associated  deflection  coefficients  for  all  stator 
blade  angle  settings.   Figure  46  is  a  summary  plot  of  maximum  pressure 
coefficients  for  all  stator  blade  angle  settings. 

On  each  of  Fig.  31  through  Fig.  36  the  3-stage  efficiency  reaches 
a  peak  value  for  each  stator  stagger  angle.  The  referred  flow  rate  at 
which  this  peak  occurs  decreases  as  the  stator  stagger  angle  increases, 

Also  plotted  on  each  of  Fig.  31  through  Fig.  34  are  the  results 

predicted  by  AXC03  for  the  three  estimated  values  of  C^  .   of  0.000, 

r  J  Dm in 

0.006  and  0.008.   The  predicted  3-stage  efficiency  reaches  a  peak  for 
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each  C^  .  •   The  referred  flow  rate  associated  with  predicted  peak 
Drain 

efficiency  decreases  with  increase  in  stator  stagger  angle.  The 
referred  flow  rate  at  which  the  predicted  peak  efficiency  occurs 
decreases  with  increase  in  Cn  .   at  a  given  stator  stagger  angle. 

AXC03  predicts  the  condition  of  surge  for  each  case  according 
to  the  criteria  of  Howell  shown  in  Fig.  29.   The  surge  condition  is 
assumed  to  exist  whenever  the  parameter  (C~6.0)/€.q  exceeds  the 
maximum  positive  value  SMAX  which  has  been  inserted  in  the  program. 
The  incidence  angle  for  the  rotor  blade  rows  becomes  larger  as  the 
through-flow  velocity  decreases  until  a  surge  condition  is  indicated. 
AXC03  predicts  surge  at  the  same  referred  flow  rate  for  each  setting 
of  stator  stagger  angle.   This  is  due  to  the  fixed  angle  of  the 
inlet  guide  vanes.   The  assumption  has  been  made  that  the  air  flow 
angle  leaving  a  blade  row  does  not  change  with  variations  in  flow 
rate  in  the  incompressible  flow  regime.   This  result  was  verified  by 
Bowen  (5).   The  inlet  guide  vanes  were  not  rotated  as  the  stator 
stagger  angle  was  varied.   At  all  stator  stagger  angle  settings  the 
first  stage  rotor  had  the  same  incidence  angle  at  particular  flow 
rates.   In  the  stage-by-stage  analysis,  "surge  in  rotor  1"  was 
indicated  at  the  same  flow  rate  in  each  run,  regardless  of  the  stator 
stagger  angle  setting  of  subsequent  stators.   The  region  of  referred 
flow  rates  smaller  than  that  at  which  surge  is  indicated  by  the 
program  is  marked  with  a  dashed  line  on  Figs.  31  through  34,  and  is 
less  than  41.5  lbm  y  R/sec  psia  for  every  case. 

When  the  value  of  the  parameter  (£""*• /p)/^©   ^s  determined  to  be 
less  than  SMI,  the  negative  value  of  the  parameter  which  has  been 
inserted  in  AXC03,  the  program  stops  the  computation  and  prints 
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"minimum  pressure  in  blade  row  xx."   The  referred  flow  rate  at  which 
computation  is  stopped  decreases  as  stator  stagger  angle  is  increased. 
At  stator  stagger  angles  of  31.8  and  lower  (Figs.  31  through  33)  the 
referred  flow  rate  at  which  computation  is  stopped  is  larger  than  the 
maximum  flow  rate  of  the  graph.   At  a  stator  stagger  angle  of  35.8 
(Fig.  34),  the  computation  stops  at  a  referred  flow  rate  of  45.5 
lbm "y  R/sec  psia.   At  stator  stagger  angles  of  39.8  and  44.3   (Figs. 
35  and  36),  the  program  did  not  establish  useful  data.   At  the  high 
stator  angles  the  referred  flow  rate  at  which  computation  was  stopped 
was  lower  than  the  referred  flow  rate  for  surge  in  the  first  stage 
rotor.   The  latter  flow  rate  remained  unchanged  at  a  fixed  inlet  guide 
vane  angle. 

This  difficulty  illustrates  some  of  the  limitations  of  program 
AXC03  when  applied  to  an  actual  machine.   The  program  assumes  that 
the  compressor  is  designed  properly.   Changes  from  stage  to  stage 
must  be  gradual.   For  example  the  annulus  area  may  change  gradually 
from  stage  to  stage  with  no  adverse  effects.   The  checks  for  minimum 
pressure  and  surge  limits  are  applied  in  each  blade  row.   Surge 
in  any  blade  row  is  interpreted  as  surge  in  the  machine,  and  similar 
limitations  are  imposed  for  minimum  pressure.   The  program  cannot 
cope  with  discontinuities  similar  to  those  that  occur  by  leaving 
the  angle  of  the  inlet  guide  vanes  unchanged.   It  is  possible  to 
imagine  a  situation  where  the  angle  setting  of  the  inlet  guide  vanes 
might  have  to  remain  unchanged;  for  instance,  if  the  inlet  guide 
vanes  must  support  the  front  bearing  of  the  rotor  shaft.   In  such  a 
design  the  first  rotor  would  have  higher  blade  efficiencies  than 
succeeding  blade  rows  at  large  flow  rates^while  at  lower  flow  rates 
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the  blade  efficiencies  of  the  first  rotor  would  be  smaller  than  for 
succeeding  blade  rows.   The  3-stage  efficiency  of  an  actual  machine 
is  an  overall  parameter  which  includes  the  effect  of  the  different 
blading  efficiencies  for  different  blade  rows.   In  an  actual  machine 
the  flow  will  adjust  itself  to  these  conditions, whereas  program  AXC03 
cannot  cope  with  these  peculiar  circumstances.   For  these  reasons  the 
program  did  not  produce  performance  data  for  the  higher  stator  stagger 
angles  of  39.8  and  44.3  . 

The  detection  of  the  actual  surge  point  during  the  tests  is 
difficult  with  the  available  instrumentation  and  must  be  based  on 
acoustic  phenomena.   Incipient  surge, probably  due  to  rotating  stall* 
was  associated  with  an  unmistakable  oscillating  change  in  sound  level 
produced  by  the  compressor.   This  effect  appeared  suddenly  even  when 
flow  rate  was  decreased  slowly.   To  eliminate  the  sound  due  to  this 
surge  condition  it  was  necessary  to  increase  the  flow  rate  by  about 
5  to  10  per  cent.   A  flow  rate  slightly  above  the  surge  point  could 
then  be  reached  by  gradual  throttling.   For  each  run  one  test  point 
was  taken  in  the  surge  region  and  one  test  point  as  close  as  possible 
to  the  surge  point.   Efficiency  at  test  points  in  the  surge  region  was 
below  0.90  for  all  stator  stagger  angles.   The  efficiency  at  the  data 
point  close  to  stall  rose  from  0.89  to  0.94  as  the  stator  stagger 
angle  was  increased  from  23.8   to  31.8  ,  and  efficiency  decreased  to 
0.89  with  further  increase  in  stator  stagger  angle  to  44.3  .   As  the 
surge  condition  is  approached  the  rate  of  decrease  of  efficiency  is 
set  by  the  relative  location  on  the  efficiency  versus  referred  flow 
rate  graph  of  the  two  data  points  in  the  surge  region  and  close  to 
surge.   At  a  stator  stagger  angle  of  31.8  which  produced  the  maximum 
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3-stage  efficiency  and  at  35.8  ,  efficiency  decreased  most  rapidly  as 
the  surge  condition  was  approached.   For  other  stator  stagger  angles 
the  efficiency  decreased  less  rapidly.   It  is  not  possible  to  make  a 
meaningful  comparison  between  measured  surge  point  data  and  the  surge 
condition  calculated  by  AXC03  because  surge  was  predicted  to  occur  at 
the  same  referred  flow  rate  for  each  case. 

It  is  possible  to  compare  measured  values  of  maximum  3-stage 
efficiency  with  the  maxima  calculated  by  program  AXC03.   For  all 
four  stagger  angles  at  which  computations  were  made,  the  maximum 
3-stage  efficiencies  were  predicted  to  occur  at  flow  rates  about 
5  per  cent  greater  than  the  flow  rates  at  which  the  measured  maxima 
occurred.   According  to  Howell  the  point  of  maximum  efficiency 
occurs  for  a  value  of  the  parameter  (t~^,0  )/£/o  about  equal  to 
+0.19  (Fig.  29).   At  greater  positive  values  of  the  parameter  the 
losses  are  believed  to  increase  rapidly.   Since  the  measured  maximum 
efficiency  occurred  at  a  flow  rate  which  was  about  5  per  cent  less 
than  was  predicted,  an  analysis  of  off -design  performance  of  a  com- 
pressor stage  was  made.   The  method  of  Vavra  (Ref.  1)  assumes  that 
the  air  flow  angles  leaving  a  blade  row  do  not  change  with  flow  rate 
in  incompressible  flow.   The  velocity  triangle  changes  caused  by  a 
reduction  from  design  flow  rate  are  shown  in  Fig.  47.   The  subscript 
d  refers  to  the  design  condition,  and  the  primed  quantities  refer  to 
off-design.   The  analysis  was  conducted  at  the  maximum  efficiency 
point  for  a  stator  stagger  angle  of  27.8  .   The  following  quantities 
were  measured; 

Design  flow  coefficient    Cp^    =  0.517 

Inlet  flow  angle  o(       (fixed)    =  20.0° 


54 


The  off-design  flow  rate  was  fixed  at  5  per  cent  less  than  the  design 
flow  rate, to  correspond  to  the  flow  rate  at  which  the  measured  maximum 
efficiencies  occurred.   The  off-design  flow  coefficient  is 
4>       =  (0.95)  (0.517) 


The  peripheral  components  of  velocity  are: 


■^     D.?5"  <t>  tan  «, 


The  relative  inlet  flow  angles  for  both  conditions  are: 

v '  0.15  <p,        L°^ff+d      J 

With  the  same  rotor  stagger  angle,  the  difference  in  incidence 
angle  I     between  design  flow  rate  and  off-design  flow  rate  is  the 
difference  in  6^ 

t'-ti  *  ?.'-  P.*  =  ***•-«*•  -  +l7° 

The  difference  between  the  parameter  ( i- C  ja  ) /£ )0   at  the  design 
flow  rate  and  the  parameter  at  a  flow  rate  5  per  cent  less  than 
design  is 

A  U-l,o)    _    ^.7    =  ^  0|5 

where  6|0  was  calculated  by  AXC03.   The  measured  maximum  efficiencies 
occurred  at  a  value  of  the  parameter  (  t- ^  iQ  )/  £    Q      which  was  greater 
by  0.15  than  the  value  of  the  parameter  for  which  the  program  AXC03 
predicted  maximum  efficiencies. 

Lieblein's  data  were  obtained  from  compilations  and  reductions  of 
cascade  test  data  primarily  for  blade  profiles  which  were  NACA 
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65  (A, n)  -  series  airfoils  built  up  on  equivalent  circular  arc 

camber  lines.   The  free- vortex  blades  installed  in  the  compressor 

have  a  thickness  distribution  which  was  derived  from  the  theory  of 

Ref.  1  applied  to  a  parabolic  mean  camber  line.   If  the  actual 

incidence  for  minimum  profile  loss  for  the  installed  blades  were 

smaller  by  1.7   than  calculated  by  subroutine  THEORYs  all  values  of 

the  parameter  (l-6.0)/£|0  calculated  in  subroutine  CASCAD  would  be 

larger  by  0.15  than  they  actually  are.   The  flow  rate  at  which  the 

maximum  efficiency  is  calculated  would  be  5  per  cent  greater  than 

the  flow  rate  at  which  maximum  efficiency  was  actually  measured. 

Howell  presented  results  of  cascade  tests  of  British  C-l  and 

C-2  airfoil  shapes  for  several  values  of  solidity.   The  C^/C„  . 

D  Dmin 

curve  of  Fig.  29  was  obtained  by  interpolation  for  a  solidity  of 

0.8.   If  the  upper  limit  of  the  parameter  ( C- C ^ 0)/€ ,0   for  minimum 

loss  were  extended  by  0.15  from  0.19  to  0.34,  the  flow  rate  at  which 

maximum  efficiency  was  calculated  analytically  would  coincide  with 

the  flow  rate  for  measured  maximum  efficiency. 

The  measured  efficiency  curves  are  more  peaked  than  the  curve 

computed  by  AXC03.   If  the  C_/C   .   curve  of  Fig.  29  had  a  greater 

slope  in  regions  where  C^/C,,  .   approaches  a  value  of  2.0,  the  pre- 

D  Dmin 

dieted  efficiency  curves  would  be  less  flat  and  would  more  nearly 

approximate  the  curves  of  measured  performance. 

It  was  one  object  of  this  thesis  to  attempt  to  determine  the 

value  of  C   .   for  the  installed  blade  profile.   Vavra  (1),  p.  377, 
Dmin 

indicates  that  these  values  usually  do  not  exceed  0.008.   For  each 

stator  stagger  angle,  analytical  calculations  were  made  for  estimated 

C  of  0.000,  0.006  and  0.008.   The  results  for  the  smaller  stator 

Dmin 
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stagger  angles  are  plotted  on  Fig.  31  through  Fig.  34.   The  value  of 

C~  .   was  estimated  for  each  stator  stagger  angle: 
Dm  In 


Stator  Stagger 

Angle 

Estimated  C^    . 

Drain 

23.8 

0.008 

27.8 

0.004 

31.8 

0.003 

35.8 

0.007 

The  arithmetic  average  of  the  estimated  values  is  0.0055,  but  there  is 
is  obviously  considerable  scatter  in  the  experimental  data. 

A  summary  plot  of  maximum  peak  efficiencies  is  presented  in  Fig. 
37.   If  continuous  control  of  stator  stagger  angle  were  possible,  as  in 
the  proposed  SST  engine,  Fig.  37  would  be  an  operating  curve  for  max- 
imizing efficiency  at  any  flow  rate-   The  peak  maximum  3-stage  efficien- 
cy was  0.956,  measured  at  a  stator  stagger  angle  of  31.8.   Also  shown  is 

a  flatter  curve  of  predicted  peak  3-stage  efficiencies  using  C     of 

Dmin 

0.006  for  the  four  lower  stator  stagger  angles.   The  maximum  predicted 
value  is  0.942  at  a  stator  stagger  angle  of  30  .   The  predicted  values 
occur  at  flow  rates  5  per  cent  higher  than  corresponding  measured  val- 
ues.  The  measured  values  decline  more  steeply  at  stagger  angles  other 
than  the  angles  for  maximum  efficiency  than  do  the  analytically  calcu- 
lated values o   The  curve  relating  stator  stagger  angle  to  performance 
shows  that, in  the  region  near  the  maximum  efficiency,  relatively  small 
changes  in  angle  are  required  for  large  variations  in  flow  rate. 
Between  stator  stagger  angle  settings  of  27.8  and  35.8  ,  a  10  per  cent 
variation  in  flow  rate  requires  only  an  8  angle  change.   In  the 
regions  more  distant  from  the  maximum  efficiency  point,  large  angle 
changes  are  required  to  handle  small  flow  rate  changes. 

Figures  31  through  36  also  contain  plots  of  3-stage  pressure  ratio 
versus  referred  flow  rate  for  each  stator  stagger  angle. 
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As  expected,  the  measured  maximum  pressure  ratio  occurred  at  a 
referred  flow  rate  between  that  for  maximum  efficiency  and  for  surge. 
It  was  not  possible  to  determine  maxima  for  the  calculated  data  of 
AXC03  because  for  each  case  the  largest  pressure  ratio  occurred  at 
the  flow  rate  tagged  "surge."  The  program  does  not  compute  perform- 
ance at  lower  flow  rates.   The  check  for  surge  occurs  in  subroutine 
CASCAD.   Based  on  cascade  data^Howell  stated  that  unacceptable  losses 
occur  at  values  of  S/£|0  greater  than  1.25.   This  conservative 
criterion  is  particularly  applicable  in  a  design  where  a  wide  range 
of  flow  rates  is  not  required.   However  for  the  purpose  of  evaluating 
pressure  ratio  for  an  existing  machine,  the  criterion  appears  to  be 
too  restrictive.   It  would  be  more  advantageous  if  AXC03  used  the 
criterion  to  predict  surge  but  then  proceeded  to  complete  the  calcu- 
lations at  that  flow  rate  rather  than  stopping  computation. 

The  measured  pressure  ratio  was  in  general  by  0.25  per  cent 
greater  than  predicted  by  AXC03,  and  the  measured  A  P     was  about 
11  per  cent  greater.   Total  pressure  is  relatively  simple  to  measure. 
The  yaw  and  pitch  errors  of  the  3-hole  probe  are  very  small.   Any 
error  in  total  pressure  measurement  is  expected  to  produce  a  total 
pressure  lower  than  expected.   Since  the  predictions  are  lower  at 
every  flow  rate,  it  is  possible  that  the  difference  in  blading 
between  the  installed  blade  profiles  and  the  NACA  airfoils,  which 
are  represented  by  Liebleins's  data,  may  be  the  cause  of  the  differ- 
ence as  was  suggested  in  the  discussion  of  3-stage  efficiency.   A 
summary  plot  of  maximum  3-stage  pressure  ratio  on  an  expanded  scale 
is  presented  in  Fig.  38.   The  peak  value  of  maximum  pressure  ratio 
occurred  at  a  stator  stagger  angle  of  31.8 
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While  writing  the  thesis,  an  error  was  discovered  in  Eq.  1,  which 

converts  the  reading  of  the  mercury  barometer  P    in  in.  of  mercury 

Bar 

to   P4    in   psfa  and  corrects    for   the    temperature  variation  of   the 
A 

specific  gravity  of   the   liquid.      The   incorrect  equation 

PA  =     PB^.   (0.492.&   +  TBir     0.5/lOOOo)l44      (psfa) 

was  used.      At   70.0  F   the  constant  by  which  P  was  multiplied  was 

Bar 

71.457  psfa/in.    of  Hg.      Reference   10  lists   the   specific  gravity  of 

"      "  00 

mercury  sy  for   the   temperature  range  0  F  to   150  F  as 

sy  =      /o.fc38  -    1.  35"4x  l0"S(TBar)  (g/cm3) 

The  conversion   factor  is 

Ci    =        (ol.&^Z  5y  .  (psfa/in.    of  Hg) 

At  70  F  the  correct  conversion  factor  is  69.55,  or  the  atmospheric 
pressure  on  which  all  calculations  were  based  was  3  per  cent  higher 
than  the  actual  atmospheric  pressure.   By  Eqs.  3,  4  and  7  the  cal- 
culated static  pressure  P  and  air  density  O      at  the  permanent 
Pitot-static  tube, and  the  total  pressure  after  the  inlet  guide  vanes 

P  .   twere  3  per  cent  higher  than  actual  also.   By  Eqs.  6,  8  and  10 
tiga 

the  calculated  velocity  V  ,  volume  flow  rate  \p c  ,  and  referred 

flow  rate  w  were  0.983  of  the  actual  values.   The  calculated  pres- 
r 

sure  ratio  P     was  0.998  of  the  actual  value  by  Eq.  11.   By  Eqs.  9 
and  14  the  values  of  weight  flow  rate  w  and  3-stage  efficiency 
were  calculated  to  be  1.2  per  cent  higher  than  the  actual  values. 

The  test  results  are  presented  in  non-dimensional  form  in 
Figs.  39  through  44.   Plotted  against  the  flow  coefficient  <^>Av ,  the 
efficiencies  follow  the  same  pattern  as  when  plotted  against  referred 
flow  rate, since  Eqs.  10  and  13  differ  only  by  a  constant.   Figure  45 
is  a  summary  plot  of  the  maximum  3-stage  efficiciency  for  each  stator 
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stagger  angle  versus  flow  coefficient  and  is  similar  to  Fig.  37  for 
the  same  reason.   Figure  46  is  a  summary  plot  of  maximum  pressure  rise 
coefficient  for  each  stator  stagger  angle.   The  pressure  rise  reaches 
a  maximum  at  a  blade  angle  of  33  and  a  flow  coefficient  of  0.475.   The 
curve  of  corresponding  stator  stagger  angles  has  an  almost  constant  slope. 

Theoretically  the  deflection  coefficient  Tw  represents  the 
change  in  peripheral  flow  components  in  a  blade  row  divided  by  the 
peripheral  speed.   Figure  47  is  a  stage  velocity  triangle  showing 
design  and  off-design  conditions  in  a  compressor  stage.   At  the  design 
condition 

d      =   U    =    U 

An  assumption  in  this  analysis  is  that  the  air  leaving  angle  from  a 
blade  row  does  not  vary  with  incidence.   Hence  the  angles  o<  j   and  /3 
remain  constant.   The  relation  between  the  off-design  deflection 
coefficient  L        and  the  off-design  flow  coefficient  Y       is  fixed  by 
similar  triangles 

1-6  <f> 


'  "  C.         + 


d 
and 


t  =   /  -  ( i-  rj 


<*>' 


or   ^.   is  expected  to  vary  linearly  with  <P     . 

At  flow  coefficients  greater  than  that  for  maximum  3-stage 
efficiency,  calculated  deflection  coefficients   c  As/  varied  linearly 
with  flow  coefficient  in  Figs.  39  through  44.   For  most  runs  the 
graph  of  L   ay   broke  sharply  from  the  linear  segment  at  the  maximum 
efficiency  point.   For  each  case  the  curve  continued  at  a  reduced 
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slope  in  the  region  between  the  break  and  prior  to  surge.   Experience 
seems  to  indicate  this  identifies  a  region  of  incipient  stall.   As 
the  blades  become  fully  stalled,  the  nature  of  the  flow  changes 
drastically.   It  has  been  suggested  that  this  change  is  reflected  in 
a  sharp  increase  in  the  slope  of  the  deflection  curve.   The  author 
was  reluctant  to  operate  the  compressor  in  the  region  of  stall  for 
extended  periods  of  time  because  of  related  vibration  and  noise  prob- 
lems.  However  a  slight  upturn  of  the  deflection  coefficient  curve 
was  apparent  at  stator  stagger  angles  of  23.8  and  27.8   (Figs.  39 
and  40). 

The  deflection  coefficient  X     is  a  measure  of  work  input.   It 
will  be  necessary  to  make  an  approximation  for  the  theoretical  power 

P„,  ,  •   Equation  11  can  be  rewritten 
Wth 

p 

PrVS      =        I       +  3"S  (Psfa) 

Making  a   first  order  expansion,   Eq.    13  can  be   rewritten 


(ft-lbf)/sec 


=  vvjV^    tJ  =     yr    AP  (ft-lbf)/.ec 

p  3  3  v  c  *-S 

Substituting  P ,  .  in  Eq.  12, the  efficiency  becomes 

Wth  ' 


Ivs  = 


3-5 


\%    CO 


Rewriting  Eq.    15 


61 


and  introducing  Eq.  16  into  Eq.  17,  there  is 


'*% 


H17*.*  AR  *, 


<W 


or  the  deflection  coefficient  is  a  direct  function  of  the  torque 

T   or  work  input, 
rq 

Figure  45  is  a  graph  of  the  I*  a<(/  which  was  calculated  at  max- 
imum 3-stage  efficiency  for  each  stator  stagger  angle.   The  minimum 
of  these  deflection  coefficients  occurs  at  the  maximum  efficiency. 
At  maximum  efficiency  the  minimum  work  input  was  required. 


62 


SECTION  7 
CONCLUSIONS  AND  RECOMMENDATIONS 

It  has  been  shown  that  maximum  efficiency  can  be  achieved  for 
each  flow  rate.   The  computer  program  AXC03  can  be  used  to  predict 
efficiency  within  2  per  cent^with  an  error  of  5  per  cent  in  flow 
rate.   Pressure  ratio  may  be  predicted  within  0.25  per  cent. 

It  is  recommended  that  further  measurements  be  made  at  increased 
rotor  stagger  angles.   Inlet  guide  vane  angles  should  be  varied  to 
simulate  the  effects  of  preceding  blade  rows.   Blade  shapes  with 
higher  blade  loadings  should  be  tested. 

Scatter  in  the  data  may  be  caused  by  variations  in  inlet  flow 
and  accentuated  by  oscillations  of  the  pressure  gauge.  A  proposed 
improved  inlet  duct  is  described  in  Appendix  E. 
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FIGURE    I 


COMPRESSOR     INSTALLATION 
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FIGURE    5  ASSEMBLY    OF    ROTOR    IN     LOW^K     CASiNG 
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FIGURE    12        EXTERNAL     STAGGER  ANGLE  ADJUSTMENT 
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FIGURE   13        INLET     BELLMOUTH 
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FIGURE    14        INLET     PITOT- STATIC    TRAVERSE 
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FIGURE  16         THREE-HOLE    PROBE    DETAIL 
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FIG.  17     PERMANENT      PITOT- STATIC     TUBE 
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FIGURE  18 


PERMANENT      P1T0T  -  TUBE    AND  THREE 
HOLE    PROBE     AHEAD     OF   FIRST  ROTOR 
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FIGURE   19  PRESSURE     READOUT    INSTRUMENTAT 
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FIGURE   22        TRAVERSE    CARRIAGE  AND  THREE  — 
HOLE    PROBE    BEHIND  THIRD   STATOR 
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FIGURE   24     TORQUE   METER    DURING    CALIBRATION 


FIGURE    25        TRAVERSE    OF  INLET    DUCT  RUN   D-l 
VOLUME  FLOW  RATE  512  FT3 /SEC 


75 


o  50 

(0 


u 

o 
•  25 

> 


^>0-0-CM 


Left 


15         10 


Horizontal 
Traverse 


5  0  5 

Radius,  in. 


10 


15 


Right 


75 


50 


o 

CD 


£  25 
u 

o 


Top 


Vertical 
Traverse 


15         10         5  0  5 

Radius,  in. 


10 


15 


^ 


Bottom 


88 


FIGURE  26        TRAVERSE    OF   INLET  DUCT   RUN  0-7 
VOLUME    FLOW  RATE   310  FT3 /SEC. 
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FIGURE  28     TEMPERATURE-ENTROPY   DIAGRAM  OF 
COMPRESSION    PROCESS   IN   3- STAGE 
MACHINE 
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FIGURE      30 

3- STAGE     EFFICIENCY    VERSUS    REFERRED 
FLOW    RATE  PREDICTED  BY  AXCO  3 

ROTOR  STAGGER  ANGLE  =  43.8° 
VARIOUS   STATOR   STAGGER    ANGLES 
CDmin  =0020 

W  (Ibm/sec)^  (^R")/R.  (psial 
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FIGURE  3f    3-STAGE   PRESSURE 
RATIO      8     3-STAGE      EFFICIENCY 
VERSUS      REFERRED    FLOW    RATE 

STATOR    STAGGER     ANGLE  =23.8* 
PREDICTION      BY  AXC03 

FOR     VARIOUS     CDmin 

Maximo   CDmjn  =-0-0.000 

-P-0.006 

-A-0.008 
Test  Data  — o- 

W(lbm/sec)-(T^CpR)/PA(psia) 


40 


Referred  Flow  Rate 


K 


rs 


94 


-  0.80 


0.75 


35 


1. 015 


FIGURE  32  3-STAGE    PRESSURE 
RATIO  8    3-STAGE  EFFICIENCY 
VERSUS   REFERRED   FLOW  RATE 
STATOR    STAGGER    ANGLE  =27.8° 
PREDICTION    BY  AXC0  3     FOR 
VARIOUS    CDmjn 

Maximo  CDmin  =-oo.OOO 
-D-0.006 
A0.008 

Test  Data  — O— 
W  ( Ibm/tec )  'TlV>  (J*S >/R(ptia) 
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FIGURE   33     3 -STAGE  PRESSURE 
RATIO   8     3-STAGE   EFFICIENCY 
VERSUS     REFERRED   FLOW  RATE 
STATOR    STAGGER   ANGLE* 31.8° 
PREDICTION    BY   AXC03    FOR 
VARIOUS    CDmin 

Maxima    CDmin  =0-0.000 
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FIGURE  34       3-STAGE    PRESSURE 
RATIO  8    3-STAGE  EFFICIENCY 
VERSUS  REFERRED   FLOW  RATE 
STATOR   STAGGER     ANGLE  =  35.8° 
PREDICTION  BY  AXC0  3    FOR 
VARIOUS    Comin 

Maxima  Co^  =-0-0000 
-q-0.006 
-6-0.008 

Ttst  Data— o- 
W(lbm/Mc)jT^  (>f7R)/PA  (ptia) 
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FIGURE  35    3-STAGE    PRESSURE 
RATIO  A    3-STAGE  EFFICIENCY 
VERSUS  REFERRED   FLOW    RATE 
STATOR   STAGQER  ANGLE* 39.8* 

Tost  Data— o— 
W  (lbm/toc)<J¥j0(.J^R)/PA(ptia) 


0.75 


35 


40  ,ff< 

Roforrod  Flow   Roto     W 


1 


o 


46 


98 


3- Stag*  Efficiency 
095    


3- Stag* 

Pressure 

Ratio 

1.025 


1.020 


0.80 


FIGURE  36    3-STAGE  PRESSURE 
RATIO  8   3-STAGE  EFFICIENCY 

VERSUS  REFERRED   FLOW  RATE 

STATOR   STAGGEJ*   ANGLE  =44.3° 

W(  lbm/8ecKfi>0  (J^/P^psio) 

Te»t  Data— O- 
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FIGURE  37     MAXMIUM    3- STAGE   EFFICIENCY 
VERSUS    REFERRED   FLOW  RATE 
FOR    VARIOUS   STATOR  8TA00ER  ANGLES 
PREDICTION     BY  AXCO  3 

FOR     C0rnn  *   0.006 

T«st  Data  — O- 
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FIGURE   38      MAXIMUM  3-STA6E  PRESSURE 
RATIO  VERSUS   REFERRED  FLOW  RATE  FOR 
VARIOUS    STATOR  STAQ6ER   ANGLES 

W(lbm/8tc)fl^  (^VR.  (ptia) 


20* 


30 


REFERRED   FLOW  RATE  W-k^- 


40 

1 


101 


3-Stogt  Efficiency, 
^3-S 


0.85 


0.25   


FIGURE  B9      3-STAGE    EFFICIENCY  ^ 

■  3~S 

DEFLECTION    COEFFICIENT *bv ,  AND  PRESSURE 
COEFFICIENT  y,QV  .VERSUS  FLOW  COEFFICIENT  </>ov 
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FIT3URE  40      3  -  STAGE  EFFICIENCY  7J^S 

DEFLECTION  COEFFICE NT  T^AND  PRESSURE 
COEFFICIENT  \^ov  .VERSUS   FLOW  COEFFICIENT <£av 
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FIGURE  41      3-STA6E    EFFICIENCY  T)    3_s, 

DEFLECTION    COEFFICIENT  Tw  ,  AND   PRESSURE 
COEFFICIENT \f/w  VERSUS  FLOW  COEFFICIENT^ 
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FIGURE  43      3 -STAGE  EFFICIENCY  ^-s 
DEFLECTION  COEFFICIENTT^AND         Tav 
PRESSURE  COEFFICIENT^  VERSUS    \f/ 
FLOW   COEFFICIENT  d>nu  STATOR 
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B,?3-S 

VTav 

Af0 


0.30    — 


0.20 


0.20 


0.15 


0.40 


0.45  0.50  0.55 

(k      ,     Flow  Cotfflcitnt 
r<ov 


0.60 


106 


3- Stage  Efficiency, 
V   3-S 


0.95 


0.90 


© 


FIGURE    44       3 -STAGE  EFFICIENCY  ^3.3 
DEFLECTION  COEFFICIENTT^ AND      Tav 
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FIGURE   45       MAXIMUM     3-STAGE    EFFICIENCY     77 

AND  CORRESPONDING  DEFLECTION     COEFFICIENT 
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FIGURE  46      MAXIMUM   PRESSURE  COEFFICIENT 

V'ov.mox     ANO  CORRESPONDING    STATOR 
STAGGER  ANGLES    VERSUS    FLOW  COEFFICIENT^ 
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APPENDIX  A 
CALIBRATION  OF  INSTRUMENTS 

The  measuring  devices  were  calibrated  to  determine  their  accuracy 
and  precision. 

The  Berkeley  digital  counter  is  located  near  a  Hewlett-Packard 
digital  counter  in  the  control  room   At  a  given  flow  rate,  there 
was  no  difference  between  the  two  instruments.   The  Hewlett-Packard 
device  was  recently  calibrated  with  a  cesium  wave  length  frequency 
standard  in  the  Standards  Laboratory  of  the  Electrical  Engineering 
Department  of  the  Naval  Postgraduate  School.   The  speed  measured  is 
considered  accurate  to  0.5  rpm« 

The  Baldwin -Lima -Hamilton  torque  meter  was  calibrated  statically 
as  illustrated  in  Fig.  24.   A  bar  was  attached  to  the  drive  end  of 
the  shaft  and  was  clamped  to  the  stand.   A  symmetrical  lever  bar  was 
attached  to  the  rotor  end  of  the  shaft.   The  lever  arm  is  20*0  in. 
Two  weight  pans  were  adjusted  to  1.44  lb.  each  with  lead  shot  in  a 
plastic  bag.   The  weights  themselves  were  checked  on  the  Toledo 
Precision  Scales  of  the  Cascade  Laboratory.   On  the  first  run,  dif- 
ferent readings  were  obtained  while  loading  and  unloading  the  pans. 
This  problem  was  overcome  in  subsequent  runs  by  tapping  the  stand  with 
a  mallet  prior  to  recording  the  reading.   This  simulated  the  vibration 
of  dynamic  operation  which  eliminated  any  frictional  effects.   Good 
agreement  was  obtained  during  tests  where  the  loading  was  increased 
and  decreased.   The  vossibility  of  temperature  effects  on  the  strain 
gauges  was  explored  in  a  third  test.   The  torque  meter  cannot  be 
calibrated  while  the  motor  is  running,  but  it  was  calibrated 
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immediately  after  shutdown.   The  shaft  was  only  slightly  warmer  to 
the  touch  than  ambient  temperature. 

The  data  of  the  three  calibration  runs  are  shown  in  Table  A-l. 
The  results  are  plotted  in  Fig.  A-l.   Run  A-3  revealed  no  temperature 
effect.   A  fourth  run  was  made  later  for  checking  purposes  and  did 
not  reveal  any  discrepancies.   The  calibration  tests  established  a 
constant  for  the  torque  meter  of  4.1565  foot-pounds  of  torque  per 
pound  reading  of  the  meter.   The  maximum  relative  error  among  the 
runs  is  0.6  per  cent-   The  meter  may  be  read  to  within  0.03  pounds, 
or  1  per  cent  of  the  meter  reading . 

The  newly  acquired  low  range  bourdon  tube  was  calibrated  also. 
Run  B-l  was  made  by  using  the  mercury  manometer  board  in  the  control 
room  of  the  Compressor  Laboratory.  The  transsonic  turbine  test  rig 
was  used  to  provide  the  pressure  difference  for  the  higher  pressure 
range.  The  present  compressor  tests  required  pressure  measurements 
down  to  0.6  in.  of  water.  Runs  B-3  and  B-5  were  made  at  low  pres- 
sure difference  with  a  Merriam  water  micro-manometer.  Pressure  was 
supplied  from  a  static  source. 

The  data  of  the  three  calibration  runs  are  shown  in  Table  A-2. 

The  results  are  plotted  in  Fig.  A-2„   The  calculation  of  velocities 

2 
was  set  up  using  pressure  difference  in  lb/ft  .   Measured  pressures 

are  therefore  very  small  numbers.   An  inverse  calibration  constant 

was   used.      The   bourdon   tube  has  a   constant  of  240.423  counts   per 

2 
lb/ft  .   The  maximum  relative  error  among  the  calibration  runs  was 

0.12  per  cent   The  meter  may  be  read  to  within  5  counts.   The 

meter's  servo  system  tends  to  overcorrect.   In  the  manual  null  mode 

the  meter  needle  tends  to  wander  a  similar  amount.   The  five  count 


114 


reading  error  amounts  to  1  per  cent  of  the  velocity  head  reading  in 
the  inlet  duct  during  flow  rate  calibration.   The  error  is  0.5  per 
cent  of  the  velocity  head  reading  at  the  permanent  Pi tot  tube.   The 
error  is  only  0.05  per  cent  of  the  measured  pressure  difference 
across  the  three  stages. 

The  Pitot-static  tube  used  to  traverse  the  inlet  duct  during 
flow  rate  calibration  is  a  modified  Prandtl-type  tube.  Another 
Pitot-static  tube  was  attached  to  the  traverse  shaft  and  the  tubes 
compared.   Sufficient  difference  existed  between  the  two  to  justify 
further  tests.   The  low  speed  calibration  tunnel  in  the  Cascade 
Laboratory  was  modified  to  receive  the  Pitot-static  tube  from  the 
inlet  duct. 

The  data  of  the  calibration  runs  are  s'hown  in  Table  A-3.   The 
results  are  plotted  in  Fig.  A-3.   Each  point  represents  the  average 
of  four  readings.   Run  C-3  indicated  that  the  velocity  head  of  the 
Pitot-static  tube  from  the  inlet  duct  was  at  most  1.1  per  cent  less 
than  that  of  the  calibration  tunnel.   Run  C-4  indicated  that  the 
static  pressure  below  atmospheric  was  at  most  1.8  per  cent  less  than 
that  of  the  calibration  tunnel.   These  maximum  errors  occurred  at  the 
lower  pressure  differences.   It  is  probable  that  the  effect  recorded 
is  in  part  the  reading  error  due  to  bourdon  tube  oscillation. 
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TABLE  A-l 
CALIBRATION  OF  BLH  TORQUE  METER 


RUN  A- 2 


APPLIED  TORQUE 
PT-LB 

4.167 

8.333 
12.500 
16.667 
20.833 
25.000 
29.167 
33.333 
37.500 
41.667 
45.833 
50.000 
54.167 
58.333 
62.500 
66.667 
70.833 
75.000 
79.167 
83.333 


METER  READING 
LB 

0.995 

2.005 

3.000 

3.985 

4.995 

5.985 

6.980 

7.960 

8.940 

9.970 
10.925 
11.950 
12.915 
13.930 
14.905 
15.925 
16.935 
17.945 
18.910 
19.970 


RUN  A- 3 


J 


APPLIED  TORQUE 
FT-LB 

19.792 

39.583 

59.375 

77.708 

92.708 
109.375 
126.042 


METER  READING 
LB 

4.745 

9.620 
14.540 
18.445 
22.400 
26.345 
30.285 


RUN     A- 4 


APPLIED  TORQUE 
FT-LB 


19.792 

38.958 

59.792 

76.458 

93.125 

109.792 

126.458 

130.675 


METER  READING 
LB 

4.77 

9.610 
14.545 
18.510 
22.435 
26.405 
30.370 
31.340 
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FIGURE  A-l         TORQUE      METER  CALIBRATION 


Calibration  Constant 


4  1821  FT- LB  Per  Pound 
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Average    Calibration  Constant    4.1565   FT -LB  ?9r  Pound 
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RUN 

i 
PRESSURE 
INCHES  OF 

MERCURY 

i 

10.260 
8.290 
7.445 
6.065 
5.235 
4.440 
4.035 
3.320 
2.415 
1.000 


TABLE      A- 2 
CALIBRATION  OF  0-12   in.    Hg  BOURDON  TUBE 
B-l  RUN     B-3  |1  RUN 


METER 

READING 

COUNTS 


WATER 
174,300  j  13.978 
139,450:  12.459 
125,750       11.062 


PRESSURE  METER  PRESSURE 

INCHES   OF     ;   READING      INCHES  OF 
COUNTS    ;       WATER 


102,900 
88,200 
74,800 
67,950 
57,100 
40,900 
16,300 


9.853 
6.993 
6.198 
4.992 
3.996 
3.036 
2.195 
1.602 
0.785 


17.409 

15,518 

13,784 

12,276 

8,737 

7,738 

6,230 

4,993 

3,794 

2,741 

1,997 

982 


14.519 
12.914 
11.397 
9.913 
8.471 
6.989 
5.471 
4.493 
3.470 
2.401 
1.512 
0.620 


B-5 

METER 
!  READING 
COUNTS 

1  18,218 

;  16,119 

14,222 

12,368 

10,560 

8,715 

6.826 

5,613 

4,398 

2,994 

1,885 

762 


; 
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FIGURE  A-2         BOURDON     TUBE    CALIBRATION 
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TABLE  A- 3 
CALIBRATION  OF  INLET  PITOT- STATIC  TUBE 


RUN  C-3 

TOTAL  PRESSURE 
MINUS  STATIC  PRESSURE 


PROBE 
2486.25 
2157.5 
1413.75 
1002.5 

793.75 


TUNNEL 
2508.75 
2182.5 
1421.25 
1011.25 

797.5 


RUN  C-4 

STATIC  PRESSURE 
BELOW  ATMOSPHERIC 


PROBE 
2450.0 
2190.0 
1422.5 
1087.5 

795.0 


TUNNEL 
2455.0 
2190.0 
1425.0 
1012.5 
810.0 
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FIGURE  A-3    INLET    TRAVERSE     PITOT-STATIC  CALIBRATION 
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APPENDIX  B 
FLOW  RATE  CALIBRATION  DETAILS 

Particulars  of  the  flow  rate  calibration  scheme  which  may  be  of 
special  interest  to  a  restricted  number  of  readers  are  included  in 
this  Appendix.   Section  B-l  gives  details  of  the  construction  and 
testing  of  the  variable  length  capillary  damping  devices.   Section  B-2 
gives  a  listing  of  the  program  ONRFL09  and  Section  B-3  presents  sample 
output  of  program  ONRFLO. 

B-l  Variable  Length  Capillary  Damping  Devices 

The  problem  of  oscillating  pressures  at  the  traversing  Pitot- 
static  tube  in  the  inlet  has  been  described.   Pressure  oscillations 
in  the  tubing  were  damped  with  alternating  sections  of  smaller  and 
larger  diameter  flow  passages.   Figure  B-l  is  a  picture  of  such  a 
device.   The  device  consists  of  three  lengths  of  needle  stainless 
steel  tubing,  of  0.020  in.  inside  diameter,  which  are  4  in. ,  6  in., 
and  8  in.  long.   They  are  connected  as  shown  to  polyethylene  tubing 
by  brass  "T's."   By  clamping  the  tubing  at  various  positions  it  is 
possible  to  obtain  total  capillary  lengths  of  4  in. ,  6  in. ,  8  in. ,  or 
18  in.   Figure  B-2  is  a  schematic  showing  the  operation  of  the  device. 

Two  variable  length  damping  devices  were  used,  one  each  in  the 
total  and  static  pressure  lines  from  the  traversing  Pitot-static  tube 
in  the  inlet.   Two  effects  of  these  devices  were  observed  on  velocity 
head,  the  difference  between  total  and  static  pressure.   As  the  length 
of  capillary  tube  was  increased,  the  fluctuation  of  the  pressure  gauge 
was  decreased.   But  as  the  length  of  capillary  tubing  was  increased, 


122 


the  time  to  come  to  a  steady  reading  was  also  increased.   Table  B-l 
presents  the  results  of  tests  of  various  lengths  of  capillary  tubing. 

Immediately  before  taking  the  velocity  head  reading  the  traverse 
probe  was  placed  at  a  new  radial  position.   This  task  took  about  20 
sec   to  perfornio   The  best  precision  for  a  time  to  damp  of  20  sec 
or  less  was  used.   The  combination  of  a  total  pressure  capillary 
length  of  18  in.  and  a  static  pressure  length  of  8  in.  satisfied  this 
requirement. 

Later  in  the  investigation,  when  the  simplified  pressure  selec- 
tion system  was  used,  it  became  possible  to  use  both  18  in.  lengths 
of  tubing  and  have  a  time  to  damp  to  a  steady  reading  of  less  than 
20  sec.   It  appears  that  the  12-channel  pressure  scanners  increased 
the  time  required  for  the  system  to  stabilize.   This  was  probably  due 
to  admittance  of  air  from  other  channels  during  sequencing  of  the 
scanner.   However  as  the  survey  progressed,  flow  rate  measurements 
were  made  at  lower  flow  rates.   As  the  machine  was  operated  nearer  to 
the  stalled  condition,  greater  pressure  oscillations  were  observed. 
The  meter  precision  was  no  better  than  +  15  counts,  and  sometimes  as 
much  as  +  50  counts. 

B-2   Program  ONRFLO 

The  flow  rate  calculations  described  in  Section  3  were  accom- 
plished by  program  ONRFLO  which  is  included  as  Table  B-2.   Since  the 
calculations  are  straightforward  and  have  been  described  earlier, 
they  will  not  be  repeated  here. 
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The  input  format  is  as  follows: 
Card  1:  NCASE    run  number  (arbitrary) 

(I7,3F7.2)  PBAR     reading  of  barometer  (in.) 

TBAR     temperature  at  barometer  (  F) 
TO       temperature  at  inlet  screen  (  F) 
Cards  2,  3,  and  4:      R        the  23  radii  at  which  data  were 
(8FL0.2)  taken  (1775,  17. 5 0,  . . .4. 00,  0.00, 

4.00. ..17.50,17.75)  (in.) 
Cards  5  through  50:     QI       velocity  head  in  inlet  duct  (counts) 
(4F7.0)  PSI      static  pressure  with  reference  to 

atmospheric  in  inlet  duct  (counts) 
QO       velocity  head  at  permanent  Pitot- 

static  tube  (counts) 
PSO      static  pressure  with  reference  to 
atmospheric  at  permanent  Pitot- 
static  tube  (counts) 

The  radii  of  cards  2,  3  and  4  are  repeated  internally  to  account  for 
the  cross  traverse.   Cards  5  through  27  introduce  measurements  made 
at  the  data  points  for  one  traverse,  and  data  for  the  cross  traverse 
are  introduced  on  cards  28  through  50c 

B-3  Output  of  ONRFLO 

Table  B-3  is  a  sample  of  the  output  of  ONRFLO.   It  is  for  the 
final  calibration,  run  14.   The  first  data  block  is  a  list  of  input 
data.   For  each  point  on  the  two  traverses,  QI ,  PSI,  QO,  and  PSO  are 
listed.   The  second  data  block  is  a  list  of  computed  velocities  VI, 
VO,  and  VIAD.   The  third  block  shows  how  the  "no  boundary  layer" 
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flow  rate,  VFTS  is  computed  by  summing  the  product  of  VA  and  A. 
The  fourth  block  shows  how  the  flow  rate  deficiency,  VFD,  is  computed. 
It  lists  the  inputs  to  and  outputs  from  SUBROUTINE  QTFE ,  which  per- 
forms the  integration  by  the  trapezoidal  rule.  The  final  statements 
list  the  final  computed  flow  rate  VFC  and  the  final  calibration 
constant  CCC 


125 


FIGURE     B-l      VARIABLE   LENGTH     CAPILLARY    DAMPING  DEVICE 


VARIABLE       LENGTH     CAPILLARY 
DAMPING      DEVICE 
OPERATION 


PRESSURE 
IN 


PRESSURE 
OUT 


MATERIALS 

CAPILLARY:      NEEDLE      STAINLESS     TUBE     0.020"  ID. 

VALVE:      POLYETHYLENE     TUBE      8     LAB     CLAMPS. 


CAP. 
LENGTH 
DESIRED 

VALVES 

A 

B 

C 

D 

4" 

CLOSED 

CLOSED 

OPEN 

OPEN 

6" 

OPEN 

CLOSED 

CLOSED 

OPEN 

8" 

OPEN 

OPEN 

CLOSED 

CLOSED 

18" 

CLOSED 

OPEN 

OPEN 

CLOSED 

FIG.    B-2 
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PRECISION     8     TIME    TO     DAMP    FOR 
VARIOUS     LENGTHS     OF     CAPILLARY 
TUBE     IN    TOTAL    AND    STATIC     PRESSURE 
LINES. 


CAR     LENGTH    IN. 

PRECISION, 
+    COUNTS 

TIME  TO 

DAMP  TO 

STEADY, 

SEC. 

TOTAL 

STATIC 

18 

18 

1 

45 

18 

8 

5 

20 

18 

6 

8 

10 

18 

4 

10 

5 

8 

18 

5 

45 

8 

8 

5 

35 

8 

6 

10 

20 

8 

4 

10 

5 

6 

18 

5 

30 

6 

8 

10 

15 

6 

6 

10 

10 

6 

4 

20 

10 

4 

18 

15 

30 

4 

8 

25 

15 

4 

6 

30 

10 

4 

4 

40 
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APPENDIX  C 
PERFORMANCE  MEASUREMENTS 

It  was  decided  to  take  flow  measurements  at  the  mean  streamline, 
using  the  radial  survey  holes.   The  assumption  was  made  that  measure- 
ments made  at  this  point  accurately  represented  the  flow  in  the 
entire  flow  annulus.   The  truth  of  the  assumption  is  verified  in 
Section  C-l.   Section  C-2  presents  details  concerning  the  data  re- 
duction program  ONRETA. 

C-l  Verification  of  Assumptions 

The  assumption  that  measurements  made  at  the  mean  streamline 
using  the  radial  survey  holes  accurately  represent  the  flow  in  the 
entire  flow  annulus  was  verified  in  three  ways.   The  total  pressure 
measured  at  the  mean  radius  was  shown  to  be  nearly  equal  to  the 
integrated  average  of  total  pressure  at  each  radius.   It  was  demon- 
strated that  the  total  pressure  measured  at  the  peripheral  location 
of  a  radial  survey  hole  was  nearly  equal  to  the  integrated  average  of 
total  pressure  at  each  peripheral  location  in  a  blade  channel.   The 
flow  was  shown  to  be  nearly  axisymmetric. 

Concerning  the  radial  variation  of  total  pressure,  the  existence 
of  a  boundary  layer  at  the  outer  radius  in  the  inlet  has  already  been 
shown.   It  was  further  shown  that  a  boundary  layer  existed  ahead  of 
the  first  stage  and  behind  the  last  stage.   With  the  stator  stagger 
angle  at  27.8  and  the  flow  rate  at  about  400  ft  /sec.  ,  a  radial 
traverse  was  made  at  the  radial  survey  hole  labeled  37  27' 30"  in 
survey  plane  8  of  Fig.  7.   Data  points  were  established  each  0.3  in. 


135 


from  a  radius  of  11.1  in.  to  18.0  in.   After  recording  the  barometric 

pressure  P    and  ambient  temperature  T   the  following  measurements 

were  taken: 

P  ,  -  total  pressure  with  reference  to  atmospheric  after  third 

stator  (counts) 

P  .  -  static  pressure  with  reference  to  atmospheric  after  third 
sd 

stator  (counts) 
q,   -  velocity  head  after  the  third  stator  (counts) 

The  total  pressure  at  the  inner  radius  could  not  be  measured  but  was 
assumed  to  be  the  static  pressure  at  11„1  in.  With  the  relations 
listed  on  p.  35   these  measuring  data  establish  the  following 
quantities: 

Tt     =    T0     +     451.7  <°R>  (2) 

to  ° 

Pt^=     PA     +         2.40.^3  (PSfa)  0) 


Sd 


Psda  "     *h    *      240.42  (psfa)  (3) 


With  the  assumption   that    the   total    temperature  after    the   third   stator 

T    ,    is   nearly  equal    to    the   inlet    total    temperature  T 
td  to 

o.00i»78  P»u  5.8.7   .  P^l5?2x|0^(lb  ,.c2/ft4)  (4) 

Xd      04.fclO(l44)        Jto  Tt0  ' 


%J,'       Z4C%4Z3  (PSf)  <5) 
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Vj      ■  -A     — —  (ft/sec)  (6) 


1 


The  average   total   pressure  was   calculated   from   the   integral   expression 


'Ua 


r 

/o.8 


18.0 

i*  VU  ft 


2,  it  v^  ri 


£ 


Both  integrations  were  performed  by  the  trapezoidal  rule.   Figure  C-l 
is  a  graphical  representation  of  the  integration.   The  integrated 

average  total  pressure  P  ,  was  2184.20  lb/ft  while  the  total  pressure 

2 

measured  at  the  mean  streamline  was  2184.34  lb/ft  .   The  relative 

error  between  the  two  pressures  is  0.006  per  cent.   A  second  radial 
survey  was  made  at  the  radial  survey  hole  labeled  149  57 '30"  in  the 

plane  called  S.  P.  8  on  Fig.  7.   The  integrated  average  pressure  for 

2 
this  survey  was  2183.92  lb/ft  .   The  total  pressure  measured  at  the 

2 
mean  radial  position  was  2184.32  lb/ft  .   The  relative  error  between 

the  two  pressures  was  0.018  per  cent.   Since  the  pressure  measured 
at  the  mean  radius  was  nearly  equal  to  the  integrated  average  of  pres- 
sure at  each  radial  location,  mean  streamline  values  were  used  to  mea- 
sure compressor  performance. 

The  traverse  carriage  was  used  for  peripheral  surveys.   Since  it 
would  be  necessary  to  remove  the  first  stage  rotor  to  survey  in  the  axial 
plane  ahead  of  the  first  stage,  the  investigation  was  confined  to  the 
plane  behind  the  third  stage  stator  where  measurements  could  be  made  simply. 
The  locations  of  the  radial  survey  holes  with  reference  to  the  indicated 
meridional  angle  of  the  carriage  are  shown  in  Fig.  23. 
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There  are  32  stator  blades.   Each  blade  channel  is  11.25  in  peripheral 
direction.   With  the  stator  stagger  angle  at  27.8  a  peripheral  survey 
was  made  at  the  mean  radius,  14.4  in.   The  average  reading  over  the 
entire  15  of  carriage  travel  was  9256  counts.   The  reading  taken  at 
the  radial  survey  hole  marked  37  27' 30"  on  Fig.  7  was  9246  counts,  a 
relative  error  of  0.11  per  cent.   The  average  over  the  first  11.25  of 
carriage  travel, representing  one  blade  channel,  was  9218  counts,  a  rel- 
ative error  of  0.30  per  cent.   Figure  C-2  illustrates  the  results  of 
this  mean  radius  survey.   A  second  survey  was  conducted  near  the  hub 
at  a  radius  of  11.2  in.   The  average  pressure  reading  of  the  traverse 
was  9532.47  counts.   The  pressure  measurement  at  the  radial  survey  hole 
was  9569.23  counts.   A  relative  error  of  0.36  per  cent  existed  between 
the  two  pressures.   On  the  basis  of  the  surveys  made  behind  the  third 
stage  stator  the  measurements  taken  at  the  radial  survey  taps  were 
accepted  as  accurately  representing  the  flow  in  the  entire  blade 
channel. 

An  investigation  of  the  axisymmetry  of  the  flow  was  conducted. 
The  results  are  tabulated  in  Table  C-l.   The  total  pressure  behind  the 
inlet  guide  vanes, P   ,  was  measured  at  eight  different  radial  survey 
holes  in  plane  S.  P.  2  of  Fig.  7.   The  average  reading  of  P  .   was 
578.1  counts.   The  nearest  reading  to  the  average  was  573.9  counts, 
taken  at  the  hole  labeled  33  45'.   This  hole  was  used  for  all  readings 
of  P  .  .   The  relative  error  was  0.73  per  cent.   The  average  of  read- 
ings taken  behind  the  third  stage  stator  P  ,  was  10,  821.5  counts. 

The  nearest  reading  to  the  average  was  10,  810  counts,  taken  at  the 
hole  marked  37  27' 30"  on  Fig.  7.   This  hole  was  used  for  all  readings 
of  P..-   The  relative  error  was  0.10  per  cent.   Within  the  measurement 
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capability  of  the  instrument,  the  flow  was  considered  axisymmetric 
and  the  most  representative  readings  have  been  used. 

C-2  Program  ONRETA 

The  data  reduction  described  in  Section  4  was  accomplished  by 
Program  ONRETA  (Table  C-2).   Since  the  calculations  are  straight- 
forward and  have  been  described  in  the  body  of  the  thesis,  they  will 
not  be  repeated  here 

The  input  format  is  as  follows: 


Card  1: 
(2I8,4F8.2) 


Card  2: 

(F8.1,4F8.0 

F8.1,F8.2) 


NCASE    run  number  (arbitrary) 
NRUN     number  of  data  points  per  run 
PBAR     reading  of  barometer  (in»  of  Hg) 
TBAR     temperature  of  barometer  (  F) 
RSTAG    rotor  stagger  angle  (43.8  ) 
SSTAG    stator  stagger  angle  (various   ) 
TO       temperature  at  inlet  screen  (  F) 
QO       velocity  head  at  permanent  Pitot- 

static  tube  (counts) 
PSO      static  pressure  with  reference  to 

atmospheric  at  permanent  Pitot- 

static  tube  (counts) 
PTIG     total  pressure  with  reference  to  atmo- 
spheric ahead  of  first  rotor  (counts) 

DELP     pressure  rise  across  the  3  stages 

(counts) 
TACH     reading  of  the  speed  counter  (rpm) 

TRAW     reading  of  the  torque  meter  (lb) 

Cards  3  through  (NRUN  -  1):   repeat  of  Card  2  for  remaining  data  points 
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The  first  five  entries  on  cards  2  through  (NRUN  -  1)  are  the  arithmetic 
average  of  four  readings  taken  two  minutes  apart.  The  last  two  entries 
are  the  average  of  two  readings  taken  five  minutes  apart. 

The  first  section  of  program  output  is  a  printout  of  the  input 
data.   The  second  section  is  the  calculated  data  for  the  run.   The 
third  section  is  the  calculated  nondimensional  coefficients. 

Tables  C-3  through  C-13  present  the  results  of  runs  1  through  11. 
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TABLE  C-l 
INVESTIGATION  OF  AXISYMMETRY  OF  FLOW 


TOTAL  PRESSURE  BEHIND 
INLET  GUIDE  VANES 


RADIAL  SURVEY  HOLE!  READING 


26  12 '30" 

*  37°27'30" 

93°42'30" 
o, 


149"57'30" 
AVERAGE 


*  READING  CLOSEST  TO  AVERAGE 


TOTAL  PRESSURE  BEHIND 
THIRD  STAGE  STATOR 


RADIAL  SURVEY  HOLE  ;  READING 


11°15' 


22°30' 


*  33°45' 


45v 


56°15' 


67°30' 


76°45' 


88  7' 30" 
AVERAGE 


•  564.2 

; 

!  587.2 

! 

573.9 

601.5 

i 

;  557.3 
565.9 

:  611.6 
563.2 

!  578.1 


l; 

I 

i 
I 
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APPENDIX  D 
DETAILS  OF  PREDICTION  PROGRAM 

This  Appendix  contains  details  concerning  AXC03  which  may  be  of 
interest  to  a  limited  number  of  readers.   Table  D-l  is  a  listing 
of  the  main  program  and  the  seven  subroutines.   Section  D-l  is  a 
description  of  the  method  of  inserting  data.   Table  D-2  is  a  sample 
of  program  output.   The  rewrite  of  input  data,  a  sample  of  interstage 
data, and  a  sample  summary  of  one  case  are  included.   Tables  D-3 
through  D-16  are  the  summary  printouts  of  the  cases  which  were  compared 
to  measured  performance  of  the  machine.   Table  D-17  is  a  printout  of 
program  CHECK  which  verified  the  errors  in  the  polynomials  of  sub- 
routines THEORY  and  CASCAD. 

D-l  Input  Data 

There  is  no  limit  to  the  number  of  cases  which  may  be  considered, 
since  data  are  read  in  prior  to  execution  for  each  case,  and  output 
printed  prior  to  considering  the  next  case.   The  executive  routine 
reads  input  data  and  prints  it  before  computation  begins.   The  input 
data  consists  of  thirteen  cards  per  case. 
Card  1:      NUM       case  number  (arbitrary) 
(17,         GAM       ratio  of  specific  heats 

4F7.3)       RN       referred  speed  (rpm/V  R) 
WREM      initial  referred  flow  rate 

(lbmy  R/sec  psia) 
DW        referred  flow  rate  increment 

(lbm*\/  R/sec  psia) 
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Card  2: 

LMAX 

(17, 

ZETAI 

5F7.3) 

EGV 

ETD 

EEX 

Card  3: 

DIN 

(3F7.3, 

DAV 

4F7.1, 

DOU 

3F7.3) 

AEA 

AIGV 

AGV 

ADI 

ALE 

SMAX 

SMI 

Card  4: 

AAA 

(10F7.3) 

Card  5: 

BH 

(10F7.3) 

Card  6: 

BLOC! 

(10F7.3) 

Card  7: 

WDO 

(10F7.3) 

number  of  stages  (3) 
loss  coefficient  of  inlet  duct 
(assumed  0.0) 
ZIGV      loss  coefficient  of  inlet  guide  vanes 
(assumed  0.05) 

efficiency  of  exit  guide  vanes  (assumed  0.85) 
efficiency  of  diffuser  (assumed  1.0) 
efficiency  of  exit  (assumed  1.0) 
inlet  pipe  diameter,  in. 
average  blade  diameter,  in. 

outlet  pipe  diameter,  in. 

2 
annulus  area  ahead  of  IGV,  in. 

2 
annulus  area  after  IGV,  in. 

2 
annulus  area  after  EGV,  in. 

2 
diffuser  exit  area,  in. 

inlet  guide  vane  exit  angle,  degrees 

maximum  limit,  (£-6 /0)/€ f      ,  set  by  Fig-  29 

minimum  limit,  (^_t/0)/€/0  »  set  by  Fi8-  29 

an  array  of  annulus  areas  after  rotor  one, 

2 
stator  one , . . . .stator  LMAXS  in. 

an  array  of  blade  heights  of  rotor  one... 

stator  LMAXS  in. 

an  array  of  blockage  factors  for  rotor  one... 

stator  LMAX  (assumed  1.0) 

an  array  of  work  done  factors  for  stage  one.. 

stage  LMAX  (assumed  1.0) 
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Card  8: 

BLDR 

(Rotor 

CAMR 

data) 

STAGR 

(9F7.3) 

SHR 

THR 

CHR 

DIAR 

DER 

CDR 

Card  9: 

BLDS 

(Stator 

CAMS 

data) 

STAGS 

(9F7.3) 

SHS 

THS 

CHS 

DIAS 

DES 

CDS 

Cards  10  and  12  are 
Cards  11  and  13  are 


number  of  blades  In  the  row 

camber  angle,  degrees 

stagger  angle,  degrees 

shape  correction  factor,  (1.0) 

thickness  to  chord  ratio 

blade  chord  length,  in. 

mean  blade  diameter,  in. 

tip  clearance,  in 

minimum  profile  drag  coefficient  (assumed 

various) 

number  of  blades  in  the  row 

camber  angle,  degrees 

stagger  angle,  degrees 

shape  correction  factor  (1.0) 

thickness  to  chord  ratio 

blade  chord  length,  in. 

blade  mean  diameter,  in. 

tip  clearance,  in. 

minimum  profile  drag  coefficient  (assumed 

various) 
duplicates  of  8 
duplicates  of  9 
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COMPARISON  OF  AlTI./VL  VALUtS  IN  NASA  SP-36  WITH  VALUES  CALCULATED  BY  AXCOM2 

SOLIDITY  =  0.8 


ANGLt,  CtoKEES 


ZLRC  CAMBER  INCIDENCE  VERSUS  ANGLE 
ACTGAL  INCIDENCE   CALCULATED  INCIDENCE   ACTUAL 


CALCULATED 


0. 

.0 

0. 

,0 

5. 

,0000000000 

0, 

,3299999833 

10, 

,0000000000 

0, 

,6799999475 

lb, 

.0000000000 

0. 

,9899999499 

20, 

,0000000000 

1. 

,  3094994o59 

-  b, 

•OUOUOCOOOO 

1. 

,6199998856 

30, 

.0000000000 

1, 

,9299993515 

J-3, 

,0000000000 

2, 

,2299995422 

<iO, 

•OoOOOOOOOO 

2. 

,5499992371 

m5, 

,0000000000 

2, 

,6299999237 

50, 

.OOOOOCOJOO 

3, 

,  1 I99996a56 

5  *» , 

.0000000000 

3. 

,3899993896 

Lly  , 

,0000000000 

3, 

,6l9999b856 

S'J, 

,0000000000 

3. 

,6299999237 

10. 

,0000000000 

4, 

,0099992752 

-0.0009329878 
0.3080989718 
0.6169836521 
0.9257200956 
1.2343072891 
1.5427494049 
1.8510446548 
2.1591939926 
2.4671907425 
2.7750444412 
3.0827503204 
3.3903064728 
3.6977157593 
4.0049800873 
4.3120946884 


0.0009329878 

0.0219010115 

0.0630162954 

0.0642798543 

0.0756921768 

0.0772504807 

0.0789546967 

0.0708055496 

0.0828084946 

0.0549554825 

0.0372495651 

-0.0003070831 

-0.0777158737 

-0.1749801636 

-0.3020954132 


PERCENT  ERROR 

0.0932987332 

6.6366682053 

9.2671022415 

6.4929141998 

5.7780303955 

4.7685480118 

4.0909175873 

3.1751365662 

3.2473926544 

1.9418897629 

1.1938962936 

-0.0090584978 

-2.1468467712 

-4.5686721802 

-7.5335502625 


SLOPE  FACTOR  N  VERSUS  ANGLE 


t . 

otGKtEi 

ACToAL 

FALlOk 

CALCULATED  FACTOR 

ACTUAL  -  CALCULATED 

PERCENT  ERROR 

0, 

.0 

-o. 

,0389999971 

-0.0343636572 

-0.0046363398 

11.8880500793 

5, 

,0000000000 

-0. 

,0499999470 

-0.0483931080 

-0.0016068891 

3.2137775421 

lo, 

.OCOOOOOOOO 

-0. 

,0609999970 

-0.0621372499 

0.0011372529 

-1.8643474579 

15, 

.COOOOoODOu 

-J, 

,0739999413 

-0.0760737062 

0.0020737648 

-2.8023862839 

cO, 

,0000000000 

-0, 

,0079999995 

-0.0906807184 

0.0026807189 

-3.0462713242 

c°3  , 

.OOOOOCOOOO 

-0. 

,  10  19999723 

-0. 1064355373 

0.0024355650 

-2.3418893814 

:0, 

,0000000000 

-0. 

,1209949919 

-0.  1238162518 

0.0028162599 

-2.3274869919 

i  5, 

,0000000000 

-J. 

,  140999^726 

-0.  1433003545 

0.0023003817 

-1.6314764023 

10, 

,0000000000 

-0. 

,  1619999409 

-0.1653655767 

0.0033656359 

-2.0775537491 

-1  IT  , 

.OOOOOOoOOO 

-0  , 

,  1859999895 

-0.  1904897094 

0.0044897199 

-2.4138278961 

5u, 

.OOoOOOOOOO 

-0, 

,2119999528 

-0.2191501260 

0.0071501732 

-3.3727235794 

Vj>, 

.  J^uuUOOOOO 

-o, 

,2409999467 

-0.2518250346 

0.0108250380 

-4.4917144775 

t  J, 

.  OOOOOOOOOO 

-0, 

,2  /r,9499504 

-0.2889920473 

0.0119920969 

-4.3292751312 

c->, 

.ouOOOOOOOO 

-0 . 

, 32o9999o00 

-0.331 1289549 

0.0101289749 

-3.1554431915 

70, 

,0000000000 

-0, 

, 3759999b  71 

-0.3787130713 

0.0027130842 

-0.7215649486 

Mul   ;i      Ctoftclb 


ObVIATlbN  VERSUS  ANGLE 
ACionl  ucVUTlL.,   CALCULATED  DEVIATION   ACTUAL  -  DEVIATIONf 


0.0014475845 

-0.0214747824 

-0.0293298364 

-0.0258062482 

-0.0145923495 

0.0006214976 

0.0161475539 

0.0282956958 

0.0333766937 

0.0277073383 

0.0175876617 

0.0093402863 

0.0092658997 

0.0136795044 

0.0288963318 


0 

.0 

0, 

,0 

-0, 

,0014475845 

-> 

.OOOOOOOOOO 

0, 

,0599999987 

0. 

,0814747810 

i  i 

.OOOOOOOOOO 

0, 

,  1299999952 

0. 

,  1593298316 

1  5, 

.OOOOOOOOOO 

0, 

, 2099999  7o5 

0, 

,2358062267 

•OOOOOOOOOO 

0, 

.2999999323 

0, 

,3145923018 

.  i 

.OOOOOOOOOO 

0  , 

,  34499997(>2 

0. 

,  3993784785 

1 

.OOOOOOOOOO 

0, 

,5099999905 

0, 

,4938524365 

n 

.OOOOOCOOOO 

0, 

,6249944952 

0. 

,6017042994 

1 1 

.  JOOuOOOOOG 

0, 

,  7549994905 

0. 

,7266232967 

.OCOCOOOOOO 

0, 

.  3999999762 

0, 

,8722926378 

>0 

.OOOOOOOOOO 

1 

,05494S4<i54 

1. 

,0424118042 

,  , 

.oooi-oooooo 

1 , 

,2500000000 

1. 

,2406597137 

oW 

.OOOOOOOOOO 

1 

.4794995422 

1. 

,4707336426 

(■3 

.OJJOU^<jOO<J 

1, 

, 7500000000 

1, 

,7363204956 

70 

.CCOooOOOOO 

2, 

.0o4949c  948 

2, 

,041 1033630 

PERCENT  ERROR 

0.1447584033 

-35.7912902832 

-22.5614013672 

-12.2886896133 

-4.8641147614 

0.1553744078 

3.1661853790 

4.4913787842 

4.3916702271 

3.0785923004 

1.65921*0198 

0.7472229004 

0.6260744929 

0.7816859484 

1.3959579468 


Table  D-   17.     Program  CHECK, Listing  and  Output   (cont.) 
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SLOPE  FACTOR  M  VERSUS  ANGLE 

Ai\GLt,  CEGRttS       ACTUAL  l-ACTOk      CALCULATED  FACTOR    ACTUAL  -  CALCULATED  PERCENT  ERROR 

0.0  0.3099999428  0.3104162812  -0.0004163384  -0.1343027353 

5.0000000000  0.3119999766  0.3125761747  -0.0005761981  -0.1846788526 

10.0000000UOO  0. 3149999976  0.3147908449  0.0002091527  0.0663976669 

15.0000000000  0.3169999719  0.3171648979  -0.0001649261  -0.05202714  73 

20.0000000000  0.3209999800  0.3198031187  0.0011968613  0.3728539944 

25.0000000000  0.3239999413  0.3228106499  0.0011892915  0.3670653105 

3u.0G0GC0G0O0  0.3279999495  0.3262923956  0.0017075539  0.5205957294 

3b. 0000000000  0.3329999447  0.3303530812  0.0026468635  0.7948539e53 

40.0000000000  0.3379999995  0.3350979686  0.0029020309  0.8585889935 

".5.0000000000  0.3439999819  0.3406320214  0.0033679605  0.9790582657 

50.0000000000  0.3499999642  0.3470600843  0.0029398799  0.8399657607 

55.0000000000  0.3569999933  0.3544869423  0.0025130510  0.7039358616 

60.0000000000  0.3649999499  0.3630179167  0.0019820333  0.5430226710 

t>5.000G0i.JG00  0.3750000000  0.3727575541  0.0022424459  0.5979855657 

70. 0000000000  0.3679999518  0.3838108182  0.0041891336  L. 0796728134 

SLOPE  FACTOR  DI/DD  VERSUS  ANGLE/ 

ANGLt,  LLoKLLb       ACTUAL  FACTUR      CALCULATED  FACTOR    ACTUAL  -  CALCULATED  PERCENT  ERROR 

0.0  U. 0769999623  0.0754218698  0.0015780926  2.0494716552 

b.OOOOOGOOOu  0.0789999962  0.0780347586  0.0009652376  1.2218189240 

lO.OGOOGOOuOG  0.0809999704  0.0802936554  0.0007063150  0.8719941378 

15.0000000000  0.0829999447  0.0826053023  0.0003946424  0.4754729667 

20.0000000000  0.085999965/  0.0853770375  0.0006229281  0.7243353128 

25.G0000C0000  0.0H99999738  0.0890154839  0.0009844899  1.0938777924 

30.0000000000  0.0939999819  0.0939274430  0.0000725389  0.0771690011 

J5. 0000000000  0.1019999981  0.1005203128  0.0014796853  1.4506711960 

40.0000000000  C.  1109999416  0.1092003584  0.0017995834  1.6212463379 

45.0000000000  0.1209999919  0.1203768849  0.0006231070  0.5149644613 

5G. 0000000000  0 .  1 31 9999t;95  0.1344547272  -0.0024547577  -1.8596649170 

55.00000GJ0UO  0.1499999/62  0.1518397331  -0.0018397570  -1.2265043259 

tjG.OOGGOOOGJO  0.1719999909  0.1729413271  -0.0009413362  -0.5472884774 

65. 0000000000  0.1979999542  0.1981655955  -0.000  1656413  -0.08365  72051 

/O.OOOOOOOOGC  0.2269999981  0.2279177904  -0.0009177923  -0.4043137431 


Table  D-17.  Program  CHECK, Listing  and  Output  (cont.) 
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APPENDIX  E 
PROPOSED  DESIGN  OF  AN  IMPROVED  INLET  DUCT 

The  present  inlet  duct  is  deficient  because  of  the  shape  of  its 
bellmouth,  because  of  the  small  velocities  that  occur  in  it  during 
operation  which  make  accurate  determination  of  volume  flow  rates 
impossible,  and  because  it  has  no  provision  to  eliminate  effects  of 
atmospheric  gusts  or  winds. 

The  entrance  bellmouth  should  have  more  gradual  changes  in  wall 
curvature  to  eliminate  local  separations  at  the  entrance  to  the 
cylindrical  duct.   It  was  stated  earlier  that  the  velocity  head  in 
this  duct  with  36  in.  diameter  is  about  0.6  in.  of  water,  a  value  too 
small  for  accurate  measurements.   To  avoid  large  fluctuations  of  the 
read-out  of  the  Texas  Instruments  pressure  gauge,  for  small  variation 

in  flow  rate,  the  velocity  head  should  be  increased  to  about  6  in. 

2  3 

of  water  or  33  lb/ft  ,  at  a  volume  flow  rate  of  360  ft  /sec,  giving 

a  response  at  the  instrument  of  about  8000  counts.   For  an  assumed 

2 
incompressible  flow  with  a  mass  density  Q  •   of  0.002378  lb  sec  / 

4  2 

ft  ,  and  a  velocity  head  q   of  33  lb/ft  ,  the  velocity  V.  in  the 

measuring  plane  of  the  inlet  duct  would  have  to  be 


2*    %l[ 
V;   =  \1    n~.  (ft/sec)  (6) 

2 
Hence  the  required  flow  area  is  360/166  =2.16  ft  ,  corresponding 

to  a  circular  area  with  19.9  in.  diameter.   Because  of  the  small 

pressure  rise  of  the  compressor  it  is  necessary  to  convert  the 

velocity  head  at  the  measuring  plane  into  static  pressure  rise  by 

arranging  an  optimum  diffuser  between  this  section  and  the  inlet 
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pipe  to  the  compressor  that  has  a  diameter  of  36  in.   Test  data  of 
Reference  11  show  that  the  included  diffuser  angle  should  be  8  for 
the  present  application.   Larger  diffuser  angles  produce  separation 
with  associated  flow  instabilities,  and  smaller  angles  require  ducts 
with  increased  lengths  where  the  boundary  layer  growth  is  excessive 
and  reduces  pressure  recovery.   For  a  diffuser  angle  of  S     the  length 
L  of  the  diffuser  is 


Z    fen  4° 


Figure  E-l  shows  the  design  of  the  proposed  inlet  duct.   A  well 
designed  bellmouth  from  a  diameter  of  48  in.  to  the  throat  of  19.9 
in.  requires  an  axial  length  of  about  26  in. ,  giving  a  total  length 
of  the  duct  of  116  +  26  =  142  in.   Ahead  of  the  bellmouth  inlet  a 
structure  will  be  attached  to  support  several  layers  of  fine  mesh 
screen  to  eliminate  flow  disturbances  by  gusts  in  the  surrounding 
atmosphere.   The  whole  duct  is  mounted  on  a  trolley  which  can  be 
rolled  onto  the  apron  outside  the  test  cell  to  permit  simple  attach- 
ment to  and  removal  from  the  presently  installed  cylindrical  inlet 
duct.   It  is  recommended  also  that  an  additional  honeycomb  flow 
straightener  be  installed  at  the  entrance  of  this  duct. 

The  bellmouth  and  the  diffuser  could  be  molded  out  of  plastics, 
reinforced  by  fiberglass,  by  using  a  wooden  template  which  is  split 
at  the  smallest  diameter. 
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